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Abstract 

Supersymmetric extensions of the Standard Model open up the possibility for new types of CP-violation. 
We consider the case of the Next-to-Minimal Supersymmetric Standard Model where, beyond the phases from 
the soft lagrangian, CP-violation could enter the Higgs sector directly at tree-level through complex parameters 
in the superpotential. We develop a series of Fortran subroutines, cast within the public tool NMSSMTools and 
allowing for a phenomenological analysis of the CP-violating NMSSM. This new tool performs the computation 
of the masses and couplings of the various new physics states in this model: leading corrections to the sparticle 
masses are included; the precision for the Higgs masses and couplings reaches the full one-loop and leading 
two-loop order. The two-body Higgs and top decays are also attended. We use the public tools HiggsBounds 
and HiggsSignals to test the Higgs sector. Additional subroutines check the viability of the sparticle spectrum 
in view of LEP-limits and constrain the phases of the model via a confrontation to the experimentally measured 
Electric Dipole Moments. These tools will be made publicly available in the near future. In this paper, we 
detail the workings of our code and illustrate its use via a comparison with existing results. We also consider 
some consequences of CP-violation for the NMSSM Higgs sector. 


1 Introduction 


After the discovery of a signal at a mass of about 125 GeV in the LHC Higgs searches the question of the 
identification of the associated state(s) and the underlying physics remains open. While the general properties are 
consistent so far with those expected for the Higgs boson of the Standard Model (SM), a wide range of alternatives 
could equally well fit the experimental data. In particular, softly-broken supersymmetric (SUSY) extensions of the 
SM count among the appealing options to solve the Hierarchy Problem and allow for a smooth transition to 
higher energy physics (e.g. Grand Unification, neutrino physics or weakly-coupled dark matter). The scarceness of 
evidence for new physics effects in precision physics or direct searches should also be weighed by the considerations 
that SUSY-inspired models offer a SM-like decoupling regime, but also that complex mechanisms in the Higgs or 
SUSY sectors - see e.g. - may account for this relative invisibility thus far. 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM), a singlet-extension of the simplest viable 
SUSY-inspired extension of the SM |^, has raised renewed interest ever since the Higgs discovery, notably due to 
its properties in the Higgs sector, e.g. allowing for an uplift of the mass of the SM-like Higgs related to F-terms or 
to the mixing of this state with a lighter singlet . The original motivation for this singlet extension rests with the 
‘/r-problem’ of the MSSM [^, which can be solved elegantly if this /i-term is generated dynamically, via a singlet 
vacuum expectation value (v.e.v.) |^. Correspondingly, the Za-conserving version of the NMSSM - allowing only 
cubic terms in the superpotential ~ is the most studied form of this model, while more general singlet couplings 


can be justified by higher-energy considerations - see e.g. 10 11 . Another usual feature in SUSY extensions of 


the SM is R-parity, which both constrains the possibility of baryon-number violation and provides a stable SUSY 
particle, hence a dark-matter candidate. 

A troubling fact rests with the observation that several NMSSM parameters - especially in the Higgs sector 
- can take complex values, hence lead to CP-violation beyond that in the quark sector. On the one hand, 
CP-violation is known as a cosmological necessity for baryogenesis. On the other, it receives severe limits at the 
phenomenological level, from the non-observation of Electric Dipole Moments (EDM’s; see e.g. [^). In this paper, 
we aim at presenting a tool which allows to study the NMSSM with complex parameters within the framework of 
the public code NMSSMTools In a first step, we will focus on the Za-conserving version although we plan on 
a generalization to tadpole and quadratic couplings of the singlet in the future. 
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The current version of NMSSMTools allows to perform several operations in connection with the spectrum of 
the CP-conserving NMSSM: in particular, it computes radiative corrections to the Higgs and SUSY spectrum, 
calculates the widths of Higgs decays or confronts the NMSSM parameter space to theoretical - e.g. vacuum 
stability - or phenomenological - e.g. Higgs searches, H-physics - limits. Several other tools aiming at the 
calculation of radiative corrections to the NMSSM spectrum have been developed in the past few years, e.g. 
NMSSMCALC or SoftSUSY [^. The latter focuses on the CP-conserving NMSSM, while NMSSMCALC allows for 
CP-violation but specializes in corrections to the Higgs spectrum. Multi-purpose programs can also be used in 
connection to the (complex) NMSSM and allow for a number of similar manipulations, provided the implementation 
of a model-file as input: this applies to SPHEND or FlexibleSUSY Ezl - which are usually coupled to SARAH 


to produce their input file. The Higgs-sector of the complex NMSSM and its phenomenology have been previously 
considered in several earlier works 19 . 

Our goal consists in generalizing NMSSMTools to the CP-violating case. While this task remains far from 
complete, the tools which we present here already allow for numerous operations: radiative corrections to the 
SUSY and the Higgs masses are implemented - so far, only the leading double-log corrections (beyond the full 
one-loop) are taken into account at two-loop order in the Higgs spectrum - ; Higgs and top two-body decays 
are computed; phenomenological limits from LEP SUSY searches or Higgs physics are tested - the latter via an 
interface with the public tools HiggsBounds 20 and HiggsSignals 21 finally we designed a subroutine to 


estimate the EDM’s. All these routines should become available on the NMSSMTools website 13 in the near future. 


This paper is intended to serve as a presentation of the calculations implemented in our tool, as well as a short 
illustration of its uses. In the following section, we will detail the characteristics of the model under study, the 
underlying assumptions and the tree-level spectrum. The third section will present the chain of subroutines that 
we designed and the operations which they carry out. Finally, we will consider phenomenological consequences 
and compare some of our results to the predictions of existing tools, before we conclude. 


2 Model, Phase-counting and Tree-level 

In this section, we present the details of the model under consideration, our notations as well as the spectrum at 
tree-level. 

2.1 The CP-violating NMSSM 

The NMSSM is a supersymmetry-inspired extension of the SM with soft SUSY-breaking terms. It differs from the 
minimal supersymmetric extension of the SM, the MSSM, in that it includes, in addition to the two Higgs SU{2)l- 
doublet superfields Hu and Hd with opposite hypercharge ±1, a supplemental gauge-singlet chiral superfield S. 
While the couplings of this singlet may take a more complex form in the general case, we will be considering 
only the R-parity and Z 3 -conserving NMSSM here, which is characterized by the following superpotential and 
SUSY-breaking terms: 


W = Xe^^^SHu ■ Hd + + Hd ■ QL[Yd]D]i + Hd ■ 


- - M 2 e*'^" 2 u)a'u}a - + h.c. 






\Hd\^ 


( 1 ) 

WaWa- + tl.C. (2) 

ml |5p + QL[ml]Ql + 


+ \AxE*^-SHu ■ Hd + - Hu ■ Ql[YuAuWr + Hd ■ QLlYdAd^ + Hd ■ Ll[Y,A,]E^j, + h.c. 

The ‘matter’ (super)field^ Ql, U^, D‘j^, Ll, E'^ should be understood as summed over generations and the 
parameters within brackets should correspondingly be seen as (complex) matrices. denotes the usual SU{2)l 
product, b, Wa and ga stand for the U{1)y, SU{2)l and SU{3)c gauginos, respectively. In the following g', g 
and gs will denote the corresponding gauge couplings and as = ffl/dTr. While the Z 3 -conserving NMSSM offers 
the simplest solution to the /i-problem of the MSSM, the inclusion of Z 3 -violating terms can be justified from 


higher-energy considerations 10,11 and turns up as a phenomenological necessity in view of the domain-wall 


problem. Our restriction to the Z 3 -conserving lagrangian follows considerations of simplicity and our work shall 
be extended to the Z 3 -violating case in the near future: we discuss in appendix]^ how this can be easily achieved. 

The minimization of the scalar potential will generate Higgs vacuum expectation values (v.e.v.’s) so that we 
may write the Higgs (super)fields in terms of their (real and positive) v.e.v.’s s, Vd, and their charged and 


'^We will omit the “ distinguishing the superfields from their scalar component, from now on. 
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neutral components: 


S = e'-*’ f s + 
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Hu = e 




Hi 


. 7,, I K+< \ 

Hd = j (3) 


The three ‘dynamical’ phases (/j^, (/>„ and (/>d add to the ‘static’ phases appearing in the lagrangian density (Eqs(g|). 
From now on, we will make the following replacements in our notations (which amounts to a redefinition of the 
Higgs fields): 






72 


Hu 


H+ 


. Vu + 


hl+ial 

V2 , 


Hd -(r- 


Vd 




H- 


y/2 


0A = 0A + 0s + 0M + 0d 

[10] ^ [Y„]e*‘^“ 

[Y„A„] ^ [YuAu^'l^- 

0K = 0K + 30s 

[10] ^ [Yd]e*'^‘' 

[YdAd] ^ [YdAdV^-^ 

4’Ax Pi = 4>Ax + 0s + 0li + 0d 

[Ye] ^ [Ye]e*‘^^ 

[YeAe] ^ [YeAe]e*'^‘' 


V ^2 = </>A„ + 30 s 


( 4 ) 


The Yukawa matrices may be written in terms of (real and positive) matrices Yu, Yd, Yg, diagonal in flavour 
space, using unitary transformations: 


[Yu] = XlYuXl ■ [Y,] = XtYdXi 

Redefining the quark and lepton (super)fields accordingly. 


[Ye] = XlYgX<k 


Ql 


[UlXI^ 


^ Xl^Vk ; 




Ll ^ LlXI^ 




( 5 ) 


( 6 ) 


and introducing the Cabibbo-Kobayashi-Maskawa (CKM) matrix Vckm = X'^"^X'l, the superpotential of Eq§ 
now reads: 


lY = XE^^SHu ■ Hd + - Hu ■ 




YuU^n + Hd 


(UlVckm\ 
\ Dl ) 


YdD-j, + Hd ■ L^YgE^ 


( 7 ) 


Finally, we make the following assumptions to ensure minimal flavour violation in the sfermion sector: 

• ~ X‘l^[mQ]X’[ = tuq, where mg is a diagonal (and, without loss of generality, real) matrix 
in flavour space. The approximation ‘~’ only holds for a matrix proportional to the identity, in the strict 
sense, but is viable, considering that the CKM matrix is hierarchical. Note that we will assume degeneracy 
for the first two generations of sfermions. 

• = m^, Xj^[rn]^]X^ = m'jj, = to|, X^[m‘^^]X^ = m\ are assumed diagonal. 

• X'l''\YuAu]X'^^ = YuAuE'P^^, X'l}\YdAd]X‘^di = YdAdE'^^d and Xl'^[YeAg\X'^^ = YgAgE‘<^^- are also treated 
as diagonal in flavour-space. 


Consequently, the soft SUSY-breaking lagrangian of Eqj^ reduces to: 

- = -MiE^^^bb - M2E*^^Wg,Wg, - MsE^^^guga + h.c. 

+ rriQ (uj^UL + dI^Dl^ + mljU'^U'^ + £>£. -|- m| (^nINl + eI^El^ + m%E'^Ej 

+ \Hu\'^ + ml^\Hd\'^ + + 


( 8 ) 


' [XAxE^^SHu-Hd+'^AuE^-S'^ Yh.c. 

- YuAuE^^-Hu ■ QlU^r + YdAdE^^dHd ■ QlD% + YgAgE^^^Hd ■ LlE'^^ + h.c. 


Eqsj^andJ^ fully characterize the model that we will be considering from now on - note that the three latter 
terms of EqjTjas well as the second and fourth lines of Eqj^are still implicitly summed over fermion generations. 
All the phases have been explicited and reduce, at this level, to four phases in the Higgs sector - ^p\, ^pu, ^pi, ^2 
will see that the minimization conditions further constrain these, as could be expected from the ‘dynamical’ nature 
of some phases -, three gaugino phases - 0 Mi, 0 M 2 i 4>M-i ~i three sfermion phases per generation - ‘PAd^ Pa^ 
- and the CKM phase finally. Given that we will neglect the Yukawa couplings of the first two generations, only 
the sfermion phases of the third generation will intervene in practice. 


3 










2.2 The tree-level Higgs sector 

The Higgs potential collects terms from the soft lagrangian (Eqj^, F-terms from the superpotential (Eqj^ and 
D-terms from the gauge interactions. We obtain: 


Vh = + ml\S\^ + AHa ■ Ha + h.c] + + h.c] 

+ A^ [|5p {\Hu\^ + \Ha\^) + \Hu ■ i?d|"] + kA S*^Hu ■ Ha + h.c. 




[\H^\^ - \Ha\^]^ + (9) 


The neutral part reduces to: 


V^o = + ml\S\^ - XA^ + h.c.] + ^ ^ ^2|5|4 


A^I^Pdij; 


0|2 


m 




e<‘fi>^-v>As*^H°H° + h.c. 


[\h:\^ - \H°a\^]" (10) 


At tree level, the Higgs v.e.v.’s are assumed to minimize this potential. A consequence is the cancellation of 
first derivatives with respect to the neutral Higgs fields at the minimum, which provides us with the minimization 
conditions: 

rriH = As [AaCOS(/3i + Kscos{(px - (fn)] — - >?{s^ + Vd) - ^ I ^ 

Vu 4 

V (7^^ ~h 0^ 

wL = As [Aa cos ifi + KS cos{ipx - (/?«)] — - A^(s^ + vl) +--- {vl - vl) 

'^d ^ 

m| = A [Aa cos ifi + 2ks cos{(px — <.Pk)\ — ns[A^ cos y >2 + 2ks] — X^{v‘^ + (11) 

Aa sin tpi = -KS sin(vjA - ^k) 

A^sm.^p2 = - AaSiu:^! - 2Kssm(i^A - E’k) —^ = “■lA-sm(i^A - ^k.) 

K S 


Here we see that the four phases of the Higgs sector are not independent but that, on the contrary, the minimization 
conditions relate (pi and (p 2 to ipx — ‘fin, the latter being the one and only ‘observable’ phase in the Higgs sector. 
Note that ipx and (/?« intervene independently in other parts of the spectrum however. We will make an explicit use 
of the minimization conditions of Eq 11 in the following lines, replacing to|, AAsintpi and A„sin (/?2 


by their expressions in terms of the v.e.v.’s. 

The terms of Eqj^ bilinear in the charged Higgs fields, define the 2x2 (hermitian) mass-matrix of the charged- 
Higgs states: 




{Mh±) = \ As [Aacos:^! -h Kscos{ipx - ¥>«)] “ (^A^ “ y ) 

— sin/3 cos/3\ /O 0 A /—sin/3 cos/3 
cos /3 sin /3y l^O J ^ cos (3 sin/3 


TO^± = 


[Aa cos + KS cos((/?A - ^ } {v't + v'd) 


tan /3 = 


Vd 


( 12 ) 


which determines the charged Goldstone boson = — sin /3 + cos /3 and the physical charged Higgs state 

= cos /3 + sin /3 . 

Similarly, the terms bilinear in the neutral Higgs fields provide the 6x6 (symmetric) mass-matrix of the neutral 
Higgs: = 5 ( osyJxdsyAi )^ with the notation ( ) meaning that fields are frozen to their v.e.v.’s. In 
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the base (/i°, /i°, /i°, a°, a°, a°), these entries read: 


= As [^acos(^i + Kscos{ipx - <Pk)] - 


Uu ^ 

■ KS cos((/?A - ^k)] + 2 ^ 4 ^ ) 


(Ad^o );^2 = [AaCOS^Pi 

V ^2 _|_ 2 

(Mho)= As [^aC0S(/3i + «;scos(v?a - </?«)] — + -— 7^'i’d 

'^d ^ 

\ 2 , 


VuVd 


'-'a " 

(Adffo)i3 = -A?;d [Aa costal + 2kscos{(Px - Lp^)] + 2A^sr;u 

(Mho)^^ = -Xvu [Aacos(/3i + 2kscos((/ja - <Pk.)] + 2X^svd 

(Ad ^0 ) 33 = KS [Af, cos + 4ks] + A A A cos (^1 


{^Ho)i4 — 0 

(Ad^o)^g = -SAKS'Ud Sin((/JA - 


(M%o)24 - 0 

(^Ho),, = 0 

(Ad^o) 2 e = - 3 XKSVuSm((fix - <Pk) 


(13) 

(Ad^o )34 = AKSVd sin((pA - <Pk') 

(Mffo)^^ = XksVu sin((^A - <Ph.) 

(M^o}^g = 4XKVuVdSin((px - (f^) 


(■^ho )44 = As [Ax cosipi + KS cos{p>x - ‘Pk)] — 

(■^Ho)45 = [AaCOS(^i + KS COs{^x “ ^k)] 

(Mho).. = As [AaC0S(/3i + ks cos{p>x - ‘Pk)] — 

(■^Ho) 4 e = [Ax cos If I - 2ks cos{fx - v^k)] 

{M]io)^^ = Xvu [Ax cosi^i - 2«:scos(vja - fn)] 

= —SKsAf. cos f 2 + X-AA cos fi + 4ks COs(:pA — Vk)\ 

As in the charged case, the neutral Goldstone boson can be singled out via a /3-angle rotation G° = — sin /3 a° + 
cos/3 a°. The remaining 5x5 symmetric block spanning the space (/i°, hP^, a° = cos /3 a° -I- sin/3 a°, a°) may be 
diagonalized via an orthogonal matrix : 

(m]jo ) = ^diag(TO|o, i = 1,..., 5)X^“ (14) 


which defines the mass eigenstates: 


= X" hi 




■xg\l 


X,fa° 


fjO 

xg a 


^^xg,hl + xg,hl + 


Xghl + Xga° + X(al 


(15) 


We will use the second notation which allows more clarity in the identification of the components. Additionally, 
we define Xg = cosPXg and Xg = sinpxg. 

Note that the positivity of the squared Higgs-masses is a stability condition of the vacuum. Remember also 

that, at 0“* order in the electroweak v.e.v.’s, one can isolate the CP-even and CP-odd sectors and diagonalize 

their doublet subspaces via rotations of angle —(3 / (3 (the singlet states are then unmixed), which disentangles 

the ‘light’ (then fully massless) ‘SM-like’ doublet states from the ‘heavy’ states with approximate squared-mass 

2 1 2 

M\ = As [AACOsy>i + Kscos{fx — fn)] g Jg (degenerate at this order with the charged state). 


2.3 The supersymmetric spectrum at tree-level 

The whole tree-level spectrum will be treated with further details in appendix Here we simply summarize, for 
the sake of notations, the basic ingredients concerning the treatment of masses and mixings of SUSY particles. 

i) Gluinos 

The gluinos are the fermionic partners of the gluons and, as such, form a color octet. Their bilinear terms originate 
in the soft lagrangian: —9 The mass states Ga, with mass M 3 (which we assume positive), 

then relate to the eigenstates of the SUSY vector superfield ga as Ga = —lei^’^aga- The phase shift then affects 
the couplings of the gluinos to coloured matter. 
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ii) Charginos 

The charginos are composed of the charged components of the electroweak gauginos and higgsinos. Their bilinear 
terms originate from both supersymmetry-conserving and violating terms and may be cast into the following form: 




T I 0 






XC + h.c. 


X^ = {-iw 




gvd 


gvu 

Xe^vxg 


= (16) 


We may diagonalize {■Xi^-+) with the help of two unitary matrices U and V: = ?7^diag(m^±, m^±)y. 

The mass eigenvalues may be assumed real and positive without any loss of generality and the mass eigenstates 
{i = 1, 2) relate to the gauge ones as: 


Xt = +Viuhu ; Xi =Uii[-lW ) +Ui2hj^ =Ui^{-lW 


in) Neutralinos 

The neutralinos are combinations of the neutral components of the electroweak gauginos and higgsinos. Their 
bilinear terms, resulting from both supersymmetry-conserving and violating terms, form a Majorana mass matrix: 


Vx« 


3 {■^x°) Xn + h.c. 


(M^o) 


Xn = {-tb, 


hlh°) 



0 

_g—y 

Y2'^^ 

0 

M2e*‘^“2 

-S-v 

toll 



0 

-Xe“^^s 



-Xe^'^^s 

0 

\ 0 

0 

-Xe^v’^Vd 

-Xe^'f^^Vu 


0 \ 
0 


-Xe^^^Vd 


-Xe^^^Vu 
2Ke^v.s / 


(M^o)^ (18) 


being symmetric, it can be diagonalized by a single unitary matrix N according to: 

{M^o) = 7V^diag(TOj^p,i = l,...,5)iV. Without loss of generality the eigenvalues m^p can be chosen real and 
positive (remember that N is complex) and the mass eigenstates relate to the gauge ones in the following fashion: 

Xi = Niii—'ib) + Ni2{—iw^) + Nish^ + Nah^ + = Nih(—ib) + iViui(~*'*h^) + Niuh^ + Nidh^ + NigJi) (19) 


iv) Sfermions 

The scalar partners of the SM fermions receive hermitian mass matrices. Due to our assumptions with respect to 
flavour violation, the three generations decouple. We keep a generic notation although only the Yukawa couplings 
of the third generation (u = t, d = 6, e = r) will be treated as non-vanishing in practice: 



(Mfj) 

(Ml) 

{Ml) 

{Ml) 


+ Yu'fjl + i (^ - 5^) - I’d) 

^ Yu Vu - svd] 

+ Y^vl + \(^^ +g^){vl- vl) 
V Yd [Ade^'^^d Vd - su„] 


g'^ +g\ 2 


K-vi) 


mfj + Yuvl Vd) ) 

Yd[Ade-^‘^^dVd- Xe^‘^-svu]\ 
ml+Yivl + S^ivl-vl) ) 


-g'"+g" 


4— (vf-n) 

Ye[Aee^^^^Vd-Xe-^^^sVu] 


Ye[Aee-^^^^Vd-Xe^‘^^sVu]\ 

< + YM+^-^ivl-vl)) 


( 20 ) 


Each mass matrix {Ml) - F = U, D, N, E - can be diagonalized via a special-unitary matrix , according to: 
{Ml) = X^^diag{ml^,mHX^. The positivity of the squared masses to|.. is a stability condition of the vacuum. 
The mass eigenstates are then defined as: Fi = Xf^F^ + Xfj^Ff^. 

This completes this short presentation of the tree-level spectrum. More details are presented in appendix 
together with the Higgs couplings. 
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3 A short walk-through the code 

In this section, we shall describe the operations which are conducted throughout our subroutines from the per¬ 
spective of the phenomenology of the CP-violating NMSSM. 

3.1 Interface with NMSSMTools 

Before coming to the actual computations of our code, let us remind the reader that we embed it within the 
NMSSMTools package. We actually use the NMHDECAY routines to define its input. In particular, we do not alter 
the running of parameters - such as the Yukawa, gauge or soft couplings -, e.g. to the average scale of the squarks 
of third generation. We simply use the corresponding quantities as calculated by NMHDECAY as our input and 
introduce the complex phases at this level. This is justified as the renormalization group equations (RGE’s) of 
the superpotential parameters leave the phases unaffected (at least up to two-loop order). A short subroutine 
init_CPV.f defines this interface and stores all the relevant quantities within commons of the code. The case 
of the parameters A\ and A„ is somewhat more subtle: given that the phases ipi and ip 2 are not free but, in 
our approach, determined by the minimization conditions of the potential (see Eq{^, we will only be using the 
quantities Ax cos (pi and cos (p 2 as degrees of freedom in practice. Therefore, we identify the NMHDECAY input 
for Ax and as ours for Ax cos 931 and A„ cos 1 ^ 2 - The one-loop RGE’s are correspondingly corrected. The wave- 
function scaling factors for the Higgs fields are also defined slightly differently from the original implementation 
in NMSSMTools, as we shall describe in section [HT] 

Given our discussion in section the following eight phases are added as new degrees of freedom: (px, Pk.i 
4>Mi , 4>M2 ! 4>M3 , VAt J , ^A^ ■ 


3.2 Supersymmetric spectrum 


The first actual operations which are carried out in connection to the CP-violating NMSSM consist in the cal¬ 
culation of the masses of the supersymmetric matter content. Similarly to the evaluation by NMSSMTools in the 
CP-conserving case, we take into account the leading radiative corrections to the masses. In the following, we list 
the new subroutines and provide relevant information concerning the calculations which are performed. 

i) mcha_CPV.f 


The purpose of this subroutine rests in diagonalizing the chargino mass matrix (Eq 16 1 according to = 


C/^diag(m^±,TO^±)lY. Similarly to the corresponding implementation within NMSSMTools for the CP-conserving 
case, the entries of the mass matrix receive one-loop radiative corrections which are calculated in the approximation 
where mass and gauge eigenstates coincide. The corresponding effects are presented in section 4.2 of 


22 - in the 


context the MSSM and still in the CP-conserving case. Small modifications appear in the CP-violating NMSSM, 
as gaugino and higgsino scalar couplings are rotated by phase factors of and Nevertheless, the 

factors of Bi functions as well as the corrections involving gauge bosons are immune to this phase shift, so that 
only the scalar interactions resulting in a Bq function - in the approximations of 22 , this reduces to the Higgs / 
higgsino loops - are affected. Another difference with respect to ref. 22 originates from the presence of singlets 


and singlinos in the higgsino self-energies. A summary of these corrections is explicited in appendix |C.1| 

The following steps are essentially identical to their counterparts in the tree-level case, which is treated into 
details in appendix B.1.4| we define two special-unitary matrices Uq and Vq diagonalizing the hermitian matrices 
and respectively. Uq ^A4^-+) Vq is then a diagonal matrix with, in general, 

non-real entries. We thus define the unitary matrices U and V via a phase-shift of Uq and Vg, where the phase of 
the lightest state is absorbed in U while that of the heavier one is absorbed in V: the resulting chargino masses 
are real and positive. 

ii) mneu_CPV. f 

The case of the neutralinos follows the same principles as that of the charginos. The tree-level gaugino and higgsino 
masses are corrected in accordance with the one-loop effects presented in appendix C.l We then diagonalize 


the complex symmetric neutralino mass matrix according to = A^^diag(m^p,z = I,...,5)A^. For that 

purpose, we consider the 10 X 10 real symmetric matrix ^ (A4;^,o)^, which can be diagonalized 

numerically by an orthogonal matrix Nq. We then extract a special unitary matrix Nq so that Nq Nq is 

diagonal. We finally absorb the remaining phases in a phase-shift of Nq, which defines the real and positive 


neutralino masses as well as the mixing matrix N. Details are provided in appendix B.1.4 
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iii) msferm_CPV.f 

We now turn to the sfermion masses. The hermitian tree-level mass matrices are diagonalized via special-unitary 
matrices ^ according to (Xip) = X'^ldiag(m|. ,mp )X^. We remind the reader that the parameters entering 
the matrices, e.g. the top and bottom Yukawa couplings, have been run to the average squark scale. The Yukawa 
couplings of the first two generation are neglected, so that the corresponding diagonalizing matrices are trivial. 
Details can be found in appendix |B.1.3| 

We then apply 0{as) corrections to the squark squared masses (consistently with what was implemented in 
the original CP-conserving treatment in NMSSMTools). Gluons, gluinos as well as the quartic sfermion D-term 
contribute to the sfermion self energy at this order. CP-phases - here, and (fAj ~ intervene in the gluino- 
sfermion couplings leading to a Bq function. A summary is proposed in appendix |C. 2 [ 

Finally, we check the positivity of the sfermion squared masses, a vacuum-stability requirement. 


iv) mgluino_CPV.f 

mgluino_CPV. f computes the gluino mass, including the 0{as) radiative corrections, which are obtained in a 

Relevant corrections include the gluon / gluino and the 


similar manner to the discussion in section 4.1 of 22 


quark / squark loops. Complex phases again enter the couplings of gluinos to squarks. Details are provided in 
appendix |C.3| 


3.3 Higgs masses and radiative corrections 

The following series of subroutines aim at computing the Higgs masses and mixing, including full one-loop and 
leading two-loop corrections. Consistently with the original approach in NMHDECAY, we will consider the effective 
Higgs potential at the average scale of the squarks of the third generation - denoted as Q where the running 
parameters are thus evaluated. 


3.3.1 Wave-function renormalization 


Momentum-dependent radiative corrections can be included in two fashions within the effective potential evalua¬ 
tion: one may reject them to the end of the calculation, as ‘pole-corrections’, or one may take them into account 
- at least partially - into the effective lagrangian as corrections to the kinetic terms. The latter choice leads 
to wave-function renormalization factors. While the two methods are formally equivalent, they lead to slightly 
divergent results at the numerical level, as we will discuss later. Following the original approach in NMHDECAY - 
presented e.g. in appendix C of 23 or appendix C of [^ -, we decide to include the leading terms - where p 
stands for the external energy-momentum of the Higgs self-energies -, originating in fermion or gauge effects, into 
the kinetic term of the effective lagrangian. Nevertheless, since we aim at a full computation at one-loop, all the 
missing momentum-dependent parts will be added as pole-corrections (see below). 

In the general case, the modified Higgs kinetic terms involve a hermitian (non-degenerate) matrix Zh{p^) as 
follows (here and below Si denotes any Higgs field; we work in momentum space and omit the factor 1/2 which 
should appear if the considered field is real): 


ZHip%s:s, 


/ , J 


( 21 ) 


The normal procedure then consists in rotating and scaling Zh via an invertible matrix Oh in order to recover the 
identity - Zh{p^) = oIjIOh then considering the ‘new’ set of fields with standard kinetic term Si = OuijSj. 

Yet, Eqs.(C.l) of 23 or (C.9-11) of show that a clever choice of the corrections included into Zh can 
make this procedure particularly simple, as Zh would turn out to be diagonal in the base of gauge-eigenstates. 
Restricting to neutral Higgs fields, one has (with 5si,Sj denoting the Kronecker symbol): 


ZH\ij=ZH^ [^Si,hl5Sj,hl+^Si,aoJSj,al\+ZHd ^ Sj ,h° + ^ +-^S + ^Si,aO^Sj,ao] (22) 


Indeed, considering the contributions of SM-fermions to Zh (A’c = 3 is the colour factor; while using the generic 
notations u, d, e, we will be considering only the third generation fermions since we neglect the Yukawa couplings 
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of the two first families), the deviations of the diagonal scaling factors from unity reacQ 

^SM ^ ^ rud, rud) + Y^Bo{p, m,, rrie)} (23) 

,5®Mf'™Zs =0 

Similarly, in the approximation where higgsinos and gauginos are simultaneously gauge and mass eigenstates {p 
denotes the doublet higgsino mass; mg, the singlino mass): 

|^Bo(p,Mi,/r) + ^So(p,M2,m) + A^Bo(p,M,wg)| 

(24) 

5'^’^Zs = ^ {A^Bo(p,Ai,/i) + K^Bo(p,TOg,TOg)} 

The last source of corrections to Z^ is the gauge sector - note that we will be working in the Feynmann gauge. 
Yet, the corresponding contributions are not diagonal in the gauge eigenbase, but rotated by an angle /3 (or —j3, 
depending on the CP-eigenvalue) in the doublet sector. Noticing however that tan /3 > 1 in practice, we may keep 
the sin^/? term in the wave-function scaling while rejecting the remaining sin/3 cos/3 and cos^terms for later 
treatment as pole-corrections. Then: 

^gauge^^^ = - lyg^BQ{p,Mw,Mw) + ^ Sq (P, ) | 

^gauga^^^ = 0 (25) 

Jgauge^^ =0 


Note that this choice in the gauge sector differs from the default treatment by NMSSMTools in the CP-conserving 
case (see Eq.(C.9-10) of §). where, moreover, pole corrections from the gauge sector are ignored. 

Before setting ZH^,Hi,s = 1-1- 


^SM farm _|_ _l_ ^gauga Zh^,Hj,,s, One is confronted to the remaining p^- 

dependence of these coefficients (via the loop functions Bq). In the ideal case, p^ would match the Higgs squared 
masses. This, however, is impractical since several mass eigenvalues are present: keeping this p^ dependence, 
hence working with p^-dependent fields and mass-matrices, and setting this implicit dependence separately to the 
corresponding Higgs squared mass after diagonalization of the mass matrix would be possible, yet problematic 
in a numerical evaluation of the mass matrices. The choice of in the CP-conserving case rested in adding an 
artificial dependence of Zh^^Hj on ln(M^/m() - Ma standing for the mass of the heavy doublet, toj approximat¬ 
ing the SM-like Higgs mass -, so as to mimic the correct logarithmic dependence after rotation by an angle —/3 
(approximating the tree-level diagonalizing rotation in the CP-even doublet sector): however an explicit rotation 
by the angle —jS shows that this purpose is missed as only the light state receives the proper logarithmic factor; 
in the case of the heavy doublet, the factor is wrong so that the result does not really improve on neglecting 
the logarithms h\{M\/m1) altogether. Therefore, we settle for the choice which consists in freezing the external 
momentum to a scale pn = 125 GeV, allowing for a good precision in the characteristics of the SM-like Higgs 
state - the most sensitive to radiative corrections. Adequate corrections when the mass is far from this scale are 
rejected to the level of pole-corrections. A final difference with comes from the implementation of the loop 
functions: we explicitly compute the full relevant Bq’s while only included the leading logarithmic terms in 
case of large mass hierarchies. 


A summary of the wave-function scaling factors is provided in appendix D.l 
Consistently, the neutral higgs fields are rescaled as: 












Ki ^ 












- ’ '"s' Py- I “U p^ - I “d fTy - 

V Zpi V Zs a/ Zh^ V Zhj, 

so that all related quantities (e.g. the mass matrices) must be rescaled accordingly. In particular the Higgs v.e.v.’s: 


(26) 






Vd{Q) = 


Vd 




\/Zhj, 

^We use the DR scheme. Loop functions are defined in appendix [ a| 


<Q) = 




(27) 
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All these operations are carried out in the initialization subroutine init_CPV.f. 

In the charged-sector, the p^-dependent terms are typically different from those appearing in the neutral case. 
However, to keep /3 as the relevant rotation angle in the charged sector together with the v.e.v. rescaling of Eq[^ 
we will use the same wave-function scaling factors Zn^^Hd- 














Vz- 


Hd 


(28) 


and restore the appropriate dependence at the level of the pole-corrections. 


3.3.2 Effective potential 

After this discussion relative to the kinetic terms, let us turn to the Higgs potential. At tree-level, it is given by 
Eq[^ Radiative corrections can be added to this picture by considering diagrams with vanishing external momenta 
or, equivalently, the Coleman-Weinberg formula for one-loop effects. In the DR scheme and the Landau gauge, 
the effective Higgs potential reads: 


Ve«.(iL) = VT{H) + SV^,,XH) 


<5Ve«.(iL) 


647r2 


Tr 




In 


Ml 



(29) 


where the trace applies to all fields d) of the model, with (7$ depending on the Lorentz properties of $ respectively 
1, 2, —2, —4, 3 for a real scalar, a complex scalar, a Majorana fermion, a Dirac fermion and a gauge boson - and 
Ml is the bilinear (‘squared mass’) matrix of the fields, where the dependence on Higgs fields has been kepl[^ 
Note that the gauge or Z 3 symmetries are still explicitly preserved by this potential (but not, in general, by its 
minimization). On the other hand, it involves terms of dimension > 5, so that expansions of the potential in the 
vicinity of its minimum will generically break the symmetries in an explicit way. 


i) Minimization conditions and corrections to the mass matrices 
The Higgs v.e.v.’s Vu, Vd, s - of Eq[^ remember that we are considering the potential at the scale Q - are 
now supposed to minimize the full potential of Eq 29 Consequently, the minimization conditions of Eq|Il| (at 


tree-level) must be modified to account for the radiative effects. This provides the so-called tadpole equation 


Sm^TT = — 


1 / dSV^ii, \ 


\dhl/V2/ 

r 2 _ 1 \ 


Smi = — 


i5(AAsin(pi) = - 


1 / \ 

2 s \dhyV2/ 

1 / dSV^, 


i5(A«, sin (^ 2 ) = 


2Xsvd \da^/'/2 

/ \ 
\dal/V2/ 


1 


2 ks^ 


1 / \ 

2\svu \ da\l^ j 

^ Y 

s \dal/V2/_ 


(30) 


Given that the parameters w|, AAsint/Si and Af^s\n{p 2 have been replaced by their tree-level values 

(Eq[TT|) in the tree-level Higgs mass matrices (Eq|T^and[T^, the shifts of Eqp0|must be included into the corrected 
mass matrices, in addition to the bilinear terms. For the charged Higgs mass matrix, this amounts to: 


5 {M‘]j±)^^ 


dHudH+ 2vudhlly/2/ 
dHXdH+ 2vdda0jV2/ 


<5(AI^±)i2 

5{Mh±)^^ 


JPSV^ I _ddV^\ 
dHudH+ 2vddall^l 
I \ 

dHJdHt 2vddhy^j 


(31) 


^In other words, one recovers the tree-level squared mass matrix (Ad|,) when replacing the Higgs fields by their v.e.v.’s in 
“^The notation {/) means that the function / of the Higgs fields is evaluated at the Higgs v.e.v.’s. 
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and for the neutral Higgs states: 


S(Mho)23 


6 

S{Mho)^^ 


1 / 


525 V,, 

1 55V,, 

(5/10/72)2 

Vu 5 / 1 O /72 

525 V,, 

1 55V„ 

(5/10/72)2 

Vd 5 / 1 O /72 

525 V,, 

. \ 

5 / 1 O/ 725 / 1 O /72 / 

525 V,, 

1 55Ver, 

(5a0/72)2 

Vu 5 / 1 O /72 

525 V,, 

1 55V,, ' 

(5a0/72)2 

Vd 5 / 1 O /72 


W A .2 \ 1 / \ 

^ ^“^12 2\dh0jV2dhyV2/ 

Wa.2 \ _ 1 / \ 

^ "“^13 2 \ahO 7725/^0/72 / 

w A .2 \ 1 / 525V,«, 1 5,5V,„. \ 

^ ^“^33 2 \(5/i0/y2)2 sdh0jV2/ 


(32) 


6{Mho)^^ 


1 / \ 
2 \daOjV2dayV2 / 


6{Mho)^q 


1 / _ \ 
2 \da0jV2da0/V2/ 




1 / \ 
2 \dayV2daOjV2/ 




1 / 52,5V, 


2\5/i2/725</72 




/ 


525 V, 


1 55V,, \ 


13 2 \ dh°JV2da°lV2 s 5a0 772 / 


Wa .2 \ 1 / 525 V, 

S{Mho)2^ - ^ 0 , /ptq 


^(■^?.c>34=d( 


23 2 \5/r0/725a0/72^ 

525 V,, 


1 55Ve«, \ 


34 2 \ 5/iV^5aO/72 sdaH^! 


6{Mho)qq - 

5(A^^o)2g ^ 2 ^ 
5 (A^^ 0)35 = 2 ^ 


1 / 525 V,, 

1 55V„ \ 

2 \(5a772)2 

55 / 1772 / 

525 V,, 

1 55V,, \ 

dhllV2dayV2 

525 V,, 

Vd da°lV2 I 
1 55V„ \ 

5/i0/725a0/72 

525 V,, 

Vd doXll\/2 j 
Vu 55V,, \ 

dhyV2da0lV2 

525Ver, 

svd do?J\f2 j 
Vu 55Veff. \ 


5/i7725a0/72 5a0/72 / 




1 7 525 V,, 2 55V,, 2n„ 55V,, \ 

2\5/r0/725a0/72 s 5a0/72 ^ s 2 5a0/72/ 


This concludes the presentation of the general formalism and we may now describe the various contributions 
to the effective potential which are computed within our code. 

ii) mhiggstree_CPV. f 

This subroutine simply defines the tree-level mass matrices at the scale Q according to Eqsj^and[^ However, 
the corrected Higgs masses are not the only information that we want to extract from the effective potential: the 
Higgs-to-Higgs couplings are also encoded within this formalism. Therefore, and for reasons that will become 
clear when we implement the various radiative contributions to the potential, we wish to match the full effective 
potential onto the following and simpler one: 


=M2|5|2 + ^ + h.c] + Vo(|5|2) 

+ {Ml + X^p\S\^)\H^l + {Mi + X%\Sl)\Hdl + + A7e*‘"“5*2) H^-Hd + h.c.] 


(33) 


Ar< 


2 

A 5 


Ii/,, 


+ y + A3|i/„|2|i/d|2 + X4\H^ ■ Hdl 


• Hdf + (A6e*‘"‘^|il„|2 + X7e^‘^^\Hal)Hu ■ + h.c. 


This is a subset of the most general singlet -I- two doublet potential which one can write up to dimension 4 term^ 
The gauge symmetry is observed. However the Z 3 -symmetry only holds up to terms quadratic in the doublet fields 

®Note however that we do not specify VodSp) further. If only terms of dimension < 4 are kept, then the only choice would be 

Vo(|SP) = k 2 |S|V 


11 








































and is explicitly broken by the terms in the last line of Eq |33| This potential is meant as an expansion of Eq |29| 
in the doublet fields and as we mentioned before, there is no reason why the Za-symmetry should hold in such an 
expansion. The characteristics of this potential are studied in appendix and matching the tree-level expression 
of Eq[^is straightforward (see appendix E.l). 

iii) mhiggsloop_sferm_CPV.f 

With this subroutine, we start adding radiative corrections to the effective potential, here those arising from SM- 
fermion and sfermion loops. These - particularly the contribution associated to the top - are known to convey 
the dominant radiative effect and lead e.g. to a substancial shift of the squared-mass of the SM-like Higgs boson. 

The corresponding one-loop effects to the neutral Higgs mass matrix are particularly easy to include in the 
Coleman-Weinberg formalism of Eqj^ since the bilinear terms provide relatively simple matrices (refer to the 


Note 


appendices B.1.1 and B.1.3). The details of the corrections are developed in the appendices D.2.1 and D.2.3 
that we recover Eqs.(C16-18) of in the CP-conserving limit. 

The situation is slightly more complex for the charged Higgs as well as for the Higgs-to-Higgs couplings: we then 
decide to expand the potential in terms of the doublet helds, up to quartic order and match the corresponding 
expansion onto the simplified potential of Eq |33[ This amounts to an expansion in ri/MgusYi where MgusY here 
stands for any sfermion mass. The sfermion contributions to the coefficients of Eq |33| are also provided in appendix 
D.2.3 Note that this alternative approach allows for a numerical cross-check with the corrections applied to the 


mass matrix of the neutral Higgs states with the method described in the previous paragraph. 

In addition to these one-loop effects, the subroutine mhiggsloop_sf erm_CPV. f also includes two-loop effects 
of ordeij^ 0 (T(\a 5 , leading to a product of large logarithms in the fermion sector: given that we are working 
at the average scale of squark masses, the squarks are assumed to give subleading contributions. On the other 
hand, effects associated to SM fermions and gauge bosons will not introduce any additional dependence on the 
new physics phases. The corresponding effects are implemented in the approximations of 24 - see also Eq.(C.19) 
of -, i.e. only the contributions to the quartic doublet parameters \u and Xd of Eq 33 are included - note that 


contributions to or Mj leave the analysis unaffected, while contributions to e.g. A^d can be absorbed in a shift 
of the tree-level term hence only drive a displacement in the parameter space. While these contributions are of 
two-loop order, they may still affect the mass of the SM-like Higgs state by several GeV, which is why we include 
them. Comparisons to more-elaborate two-loop calculations show that this approximation works well numerically 


(at the GeV level). Two-loop effects beyond this order have been considered in 25 


iv) mhiggsloop_inos_CPV.f 

The next subroutine implements the radiative effects associated to charginos and neutralinos. Sticking to the 
Coleman-Weinberg approach, we consider the 9x9 bilinear term associated with gauginos and higgsinos (refer 
to appendix B.1.4|. Due to the large rank of this matrix, we exclusively employ the method which consists in 


expanding the potential and matching it to the simplified version of Eq|33| The corresponding results are collected 
in appendix D.2.4 Note that they differ from e.g. Eq.(C.22-24) of where additional simplifying assumptions 


had been made. 

v) mhiggsloop_gaugehiggs_CPV.f 

The contributions of the electroweak gauge bosons to the Higgs potential seem easy to include in the Landau 
gauge: see appendix |D.2.2| Yet the drawbacks of the Landau gauge are felt in the Higgs sector, where one 
then has to handle massless Goldstone bosons. The associated infrared divergences are of course purely spurious 
and disappear once confronted to momentum-dependent corrections, as already noted in 26 . Still it remains a 


technical issue to manipulate with caution. Moreover, the strategy consisting in diagonalizing the field-dependent 
bilinear matrices, which we have been employing until here, becomes impractical, even in an expansion in terms 
of the doublet fields, due to the large number of parameters and operators involved in the Higgs bilinear terms. 
Instead, we decide to employ the concurrent strategy in Higgs-mass calculations, which simply consists in a 
direct diagrammatic evaluation of the Higgs self-energies and tadpoles generated by Higgs loops. Nevertheless, 
disentangling the Higgs and gauge contributions in this approach proves quite artihcial so that we will perform 
the calculation simultaneously for both types of particles appearing in the loop. 

Explicit expressions for the gauge and Higgs one-loop contributions to the Higgs self-energies and tadpoles are 
summarized e.g. in 22 or (with different conventions for the loop functions), in the context of the MSSM, 


and the NMSSM differs only in the definition of the couplings and the presence of the singlet helds, hence leads to 
a formally comparable result. We choose to work in the Feynmann gauge as it is then possible to set the external 
momentum to 0 without generating IR-divergent logarithms. Indeed, we still aim at computing, not only the 
corrections to the Higgs masses, but also to the Higgs-to-Higgs couplings. For this, we proceed in the following 


®The conventions or 0{at^i,as, are also used in the literature. 
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fashion: after the radiative corrections to the Higgs mass matrices are evaluated at zero momentum, we subtract 


the pure-gauge contribution in the Landau gauge (for which we already know the potential from appendix D.2.2). 
The remaining ‘Higgs’ contributions to the mass matrices can then be identified with those that a Za-conserving 
renormalizable potential would produce, allowing for a reconstruction of the corrections to the Za-conserving 


parameters of the potential: this procedure is described for the CP-conserving case in 28 and is straightforwardly 


generalized to the CP-violating case. Of course, we then miss contributions to the Za-violating parameters (the 
last line of Eq 33) but, as discussed in 28 , these are subleading in the leading-logarithmic approaclQ Further 
details can be found in appendix |D. 2.5 


This completes the list of radiative contributions implemented in the effective potential. 


3.3.3 Pole corrections 

The operations described in the previous lines have provided us with mass matrices for the Higgs states where 
radiative corrections from the potential (i.e. at zero external momentum) have been included. We will now 
detail how we account for momentum-dependent corrections. These calculations are conducted in the subroutine 
mhiggsloop_pole_CPV.f. 

First, let us remind the reader that the radiative effects associated with non vanishing external momentum 
have been partially encoded into the wave-function scaling factors of paragraph |3.3.1| It is necessary to rescale 
the Higgs mass-matrices in order to account for the re-scaling of the Higgs fields: 

A /3-angle rotation in the pseudoscalar and charged sector allows to rotate away the Goldstone bosons, leaving 
us with a 5 X 5 symmetric mass matrix ( Ad^o) for the neutral sector and a DU squared-mass for the charged 


Higgs my±^. is now diagonalized according to Eq 

for the neutral Higgs, and their rotation matrix X^° 
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providing us with corrected DR squared-masses 


We then apply pole corrections to the DR squared-masses in order to evaluate the pole masses: 


2 2 DR 

rtigo = rn Si 


(i + J2SZh,X, 


H °: 

ij 


ns?s?(4'^ — ~ nsPs‘>(?'^ — 0) 


(l + cos^ /3 + 5Zh^ sin^ 0) - [Uh+h- - Uh+h- = 0)] (35) 

While ideally the Higgs self energies Ilgggg{p0 and IIh+h-{p^) should be evaluated at the pole masses, we 
approximate the latter by the DR masses. The full one-loop pole-corrections are applied. Shifts of the wave- 
function scaling factors 5ZH^,Hd,S are know from Eqs 23|25 The shifts in the Higgs self energies are provided in 
appendix |D.4| 

This concludes our evaluation of the masses in the Higgs sector. We now wish to comment briefly on the 
precision achieved in this calculation. For this, it is instructive to consider the impact of the one-loop corrections 
with respect to the situation at tree-level. For mostly-doublet states, the leading effect is driven by the top- 
quark loop and, respectively to the tree-level mass mn, can be quantified as ^ ^^2 00- In Assuming that 
Q = O(TeV), this amounts to a correction at the percent level for ttih = O(TeV), but reaching a magnitude of 
~ 100% for rriH = 0(100 GeV): this accounts for the well-known sensitivity of the SM-like Higgs mass to radiative 
corrections. Contributions at the two-loop order will involve the strong coupling gs, or the top Yukawa coupling 
again, multiplying logarithms of a similar magnitude, so that the typical effect would easily amount to ^ 30% of 
the one-loop contribution. Now, considering that we have included the leading double-logarithmic effects in the 
calculation, we can estimate a reduced uncertainty from higher orders, say at the level of ~ 10% of the one-loop 
corrections. For a Higgs mass at ~ 125 GeV, this still amounts to an uncertainty of several GeV. For a state at 
niH = O(TeV), this reduces to the permil level. The latter accuracy is treacherous however, as other sources of 
uncertainty appear e.g. in the determination of the couplings or neglected electroweak corrections entering the 
definition of the Higgs v.e.v.’s. In the outcome, one should not expect a precision under 0(1%) for the masses of 
the heavy doublet states. Gorrections to singlet states are typically smaller, since the associated couplings - A, 
K - are of order < 0(0.7) and the hierarchies between Higgs bosons and higgsinos may not be as large as those 
between SM fermions and sfermions. However, when the singlets mix significantly with doublet states, they will 
correspondingly acquire part of the larger uncertainties on doublet masses. 


^Note that, while the parameters of the effective potential will then receive contributions which are valid only at leading logarithmic 
order, this is not the case for the contributions to the Higgs masses since they are obtained directly from the diagrammatic calculation. 


13 


























3.4 Couplings, decays and constraints 

After the previous subroutines are run, one has a complete set of corrected masses and rotation matrices at one’s 
disposal. The following move consists in confronting this spectrum to physical processes. 


3.4.1 Supersymmetric and Higgs couplings 

The couplings of supersymmetric particles and Higgs bosons can result in somewhat lengthy expressions. We thus 
design two subroutines, susycoup_CPV.f and higgscoup_CPV.f , in order to evaluate and store them within the 
code: 


• The couplings of charginos / neutralinos to sfermions and SM fermions are implemented according to the 
formulae of appendix |B. 3.1 and B.3.2 for the three generations (still neglecting the Yukawa couplings of the 
two first generations). 

• The trilinear couplings of the Higgs bosons to sfermions are computed after the results of appendix |B. 2. 3[ 

• The couplings of the Higgs bosons to charginos and neutralinos are also included as presented in appendix 

[ET41 

• Finally, we calculate corrected DR Higgs-to-Higgs couplings where radiative effects from sfermions, charginos, 
neutralinos and Higgs bosons are obtained from the simplified effective potential of Eq|^ relevant formulae 
are provided in appendix |E.4[ Corrections from fermionic and gauge loops are explicitly incorporated as 

The Yukawa and gauge couplings employed here have been run to the 


given in appendix D.2.1 and D.2.2 


scale corresponding to the mass of the neutral Higgs state associated with the first index in the coupling. 


Note that the rescaling of Higgs fields in Eqs 26 and 28 is also accounted for when computing the couplings of 
Higgs bosons. 


3.4.2 Higgs and top decays 

We then adapt the existing NMSSMTools subroutines decay.f and tdecay.f - respectively computing the Higgs 
and top two-body decays in the CP-conserving NMSSM - to the CP-violating case. 

The subroutine hidecay_CPV. f calculates the Higgs widths and the dominant branching ratios. The following 
decay channels are considered: 

• decays into a pair of SM fermions: —>■ , ss, cc, bb, ti] us, ub, cs, cb, th] 

• decays into (on- and off-shell) gauge bosons: —)■ WW, ZZ, 77, Z"f, gg-, 

• decays into one Higgs and one gauge boson: Sf —>■ ZSj, —>• 

• Higgs-to-Higgs decays: 5° —)■ 5°5'°, H+H~; 

• supersymmetric decays: S° -)■ X^Xk, XjXt, Fj Fk] H+ -^x'^xl-. F*F'^. 

In the subroutine tdecay_CPV.f , we compute the following top decays: t W'^h, H'^b, Ty^. As in the 
original CP-conserving version, leading QCD corrections have been taken into account. 


3.4.3 Phenomenological tests 

We finally propose several tools to confront the CP-violating NMSSM spectrum to experimental constraints. 

checkinin_CPV. f compares the value of the neutral effective potential at the electroweak symmetry-breaking 
minimum with that at other points, e.g. for vanishing v.e.v.’s. Loop effects from the SM fermions and gauge 
bosons are included explicitely in this evaluation, while other radiative effects are encoded within the approximate 
potential of Eq |33| We also vary the dynamical phases and check whether this generates a deeper minimum. Finally, 
the minimization conditions of EqfTT] and [30] are calculated explicitly, which allows e.g. to test the naturalness of 
the squared masses ^ of the potential: they should remain of the order of the SUSY-breaking scale. 

In constsusypart_CPV . f , we generalize to the CP-violating case LEP limits on superparticle searches that 
were included in NMSSMTools for the CP-conserving case: 

• test on chargino, slepton, gluino and squark masses; 
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limits onT ^ blN, T —)■ cx°, B —)■ 6%' 


0 . 


constraint on the invisible Z-width and neutralino-pair production. 


HBNMSSM_CPV.f converts our spectrum into input for HiggsBounds 20 and HiggsSignals 21 . This allows 


to test the Higgs sector in view of LEP, TeVatron and LHC results via a call to the subroutines included within 
these public tools - note that NMSSMTools, HiggsBounds and HiggsSignals must be interfaced to make use of 
this subroutine. The chosen input mode is that employing effective couplings (see the documentation in 20 ). 
Additional widths and branching ratios are taken from the results in hidecay_CPV . f and tdecay_CPV.f . In the 
following section, we will be using the current versions HiggsBounds_4 .2. 0 and HiggsSignals_l . 3. 1, which 
incorporate all the experimental results released till december 2014. The default uncertainty on the Higgs mass 
precision is set to 3 GeV and modeled as a gaussian distribution. HiggsBounds delivers a 95% confidence level 
cut on the NMSSM parameter space relative to limits from unsuccessful Higgs boson searches. HiggsSignals 
performs a x^-test of the Higgs properties of a given spectrum based on the current experimental characteristics 
of the the signals measured at ~ 125 GeV. The default setting of version 1.3.1 includes 81 test-channels based 
on the material released by the ATLAS, GMS, GDF and DO collaborations. The output, the x^-value, provides a 
measure of the compatibility of the tested Higgs spectrum with the observed signals. While the implementation 
within HiggsSignals accounts for more involved effects, such as correlations among channels or uncertainties, the 
X^ can be grossly understood as the sum of squared deviations between theoretical predictions and experimental 
measurements, normalized for each channel to the sum of squared theoretical and experimental uncertainties for 
this channel: therefore, the smaller the x^? the more compatible the Higgs spectrum proves in view of the measured 
Higgs signals. Statistical interpretations in terms of P-values are possible, yet depend on the details of the chosen 
tests or of the definition of the statitical ensembles: in our discussion, we will conhne to the thumb rule stating 
that x^-values of the order of the number of test-channels (81 here) are regarded as competitive. The x^-value in 
the SM limit (~ 78 in practice) gives another point of reference from which one may appreciate the quality of the 
ht of a particular Higgs spectrum to the observed Higgs signals. 

We also include an alternative set of tests for the Higgs sector, based on the original subroutines of NMSSMTools. 
These collect: 


• LEP_Higgs_CPV.f : LEP limits applying on neutral Higgs bosons produced in association with Z’s - e+e —>■ 
Z* —>■ S^Z - or in pairs - e+e” —>■ Z 


qo qo _ 
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TeVatron_CHiggs_CPV.f : TeVatron limits applying on a charged Higgs boson produced in top decays 30 


• bottomonium_CPV. f : test for a light mostly GP-odd Higgs intervening in bottomonium decays - based 
on 


31 


• LHC_Higgs_CPV. f : the inclusion of LHG limits on neutral or charged Higgs searches as well as the confronta¬ 
tion to the signals at ~ 125 GeV - after 32 


- are in progress. 


Note that these routines will not be used in the next section, as we will employ the currently more complete set 
of tests performed by HiggsBounds and HiggsSignals. 

Finally, we design a subroutine EDM_CPV.f to estimate the electric dipole moments of the electron, the thallium 
atom, the neutron and the mercury atom. We essentially follow the summary in 33 - in the context of the MSSM; 
see also 34 for previous works in the NMSSM. The supersymmetric one-loop effects are mediated by charginos, 


neutralinos or gluinos and sfermions. Moreover the two-loop diagrams of the Bar-Zee type - involving a fermion 
or sfermion loop connected to the quark / electron line by a Higgs and a photon propagator - are known to convey 
a sizable effect: these are particularly sensitive to the phases appearing in the Higgs sector. Other contributions, 
mediated by dimension 6 operators, are included as well. We estimate the associated uncertainties by adding 
linearly a 10% error on effects involving no coloured particles and a 30% error on contributions involving the 
coloured sector. Additional uncertainties associated to scale-running or hadronic parameters are also incorporated. 


4 A few applications 

We shall now make use of the subroutines which we have just presented and study phenomenological effects 
associated with the GP-violating NMSSM. This will be the opportunity to test our tool and compare its predictions 
with existing results. 
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4.1 CP-conserving limit 

Setting all the phases to zero, it is possible to consider the CP-conserving case: in particular this allows to study 
how our results connect to the precision calculations implemented within NMSSMTools. Given that the input is 
common, discrepancies directly give an insight on the differences of treatment and the numerical magnitude of the 
corresponding effects. 


mj, (GeV) Precision 0 (GeV) Precision I 




Figure 1: General aspect of the Higgs spectrum in the NMSSMTools procedures with precision 0 (upper left), 
precision 1 (upper right) and precision 2 (down left-hand) and in ours (dotted lines). A G [0,0.65], k = 0.45, 
tan/3 = 8, = 125 GeV, Ma = I TeV, = -288 GeV, mf g = 1 TeV, Q = 1.5 TeV, = 200 GeV, 

2Mi = M 2 = M 3/3 = 0.5 TeV, At = —2 TeV, A^t = —1.5 TeV. In the down right-hand plot, the singlet 
composition of the two lightest CP-even states is displayed for precision 0 (green solid lines), precision 2 (blue 
solid lines) and for our calculation (dotted lines). 


i) Higgs spectrum 

We shall first consider the Higgs masses. NMSSMTools provides three levels of precision in the inclusion of the 
radiative corrections to the Za-conserving Higgs sector: 


‘Precision O’: the default one - essentially following the procedure described in appendix C of [^ - confines to 
leading logarithmic order. Momentum-dependent effects are taken into account only to the extent of wave- 
function renormalization (where the implementation is slightly different from ours: remember the discussion 
in section 3.3.1) and pole-corrections associated with the SM-fermion sector. 
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• ‘Precision 1’: a full one-loop -I- leading two-loop (to order 0{Y^i^as)) implementation without momentum- 
dependent effects. 


‘Precision 2’: a full one-loop -I- leading two-loop (to order 0{Y^\as)) implementation including momentum- 


dependent effects. It follows the work of 35 


Formally, our implementation - full one-loop including momentum-dependent corrections -I- leading two-loop 
double logarithms of order 0{Y^\as,Y^i^) - should fall somewhere between these three procedures in terms of 
precision. 

We first test our results in a region of the parameter space where k = 0.45, tan/3 = 8 , = 125 GeV, 

Ma = 1 TeV, A„ = -288 GeV, mf ^ = 1 TeV, = 1.5 TeV, mi j; = 200 GeV, 2Mi = M 2 = M 3/3 = 0.5 TeV, 
At = —2 TeV, Ab^r = —1.5 TeV and we scan over A G [0,0.65]. The Higgs masses are displayed in Figj^andj^ 
the results of NMSSMTools for precision ‘0’ (greenish colors), ‘1’ (pink colors) and ‘2’ (bluish colors) are shown as 
solid lines while our calculation corresponds to the dots (yellow to red tones). We observe a significant variation of 
the masses corresponding to the mostly-singlet states while the doublet masses are grossly constant with varying 
A. A typical NMSSM effect develops when singlet masses are close to doublet masses, as significant mixing may 
appear. In particular, when the singlet state is slightly lighter than the doublet one, the mixing tends to uplift 
the mass of the mostly-doublet Higgs. This is what occurs in this example for the GP-even sector in the upper 
range of A. In Fig[^ we see that our results fit quite closely the predictions of the procedure with precision 2, 
while larger discrepancies appear with respect to precision 0, especially at large A. 

Figj^ allows for a closer comparison among Higgs masses. For the Higgs states with mass close to ~ 125 GeV 
(upper left-hand quadrant), we note a remarkable agreement between our calculation and the masses obtained 
with the precision setting 2, while the results obtained with precision 0 are about 2 GeV off: this fact should not 
make us forget that the uncertainties affecting the Higgs mass computations (also in the setting of precision 2) are 
of the order of several GeV. However, it justifies the observation that the leading two-loop effects are captured by 
the simpler inclusion of double logarithmic terms. 

Concerning the heavy mass states, we observe in Figj^- in the upper-right and lower-left hand quadrants - that 
our results are intermediary between the calculations with precision 0 and precision 2. However, we note that the 
leading difference with precision 2 originates in the implementation of the wave-function scaling factors. Indeed, if 
we set the ‘Z-factors’ to 1 and modify the pole-corrections accordingly, we observe that our result - corresponding 
to the khaki dots in the lower-right-hand corner of Figj^- then matches that with precision 2 somewhat more 
closely (at the permil level). It is quite easy to see how the discrepancy develops between these two procedures. 
For this, let us focus on the CP-odd doublet state, where we will neglect the mixing with the singlet. In the case 
where the wave-function scaling factors are set to 1 , the squared mass of this state is - schematically: the effect 
of potential and pole corrections are encoded as i 5 p°up°i» _ obtained as (all the p^-dependent terms are treated as 
pole corrections): 

m^lno z = [Aacospi -b Kscos{ip\ - p„) -b (l -b (36) 

In the approach where the wave-function scaling factors are taken into account at the level of the kinetic terms, 
the Z-factors intervene in the calculation at several steps: first, for the scaling of the v.e.v.’s, which transforms 
the tree-level mass-matrix in the CP-odd doublet sector to: 


As 

7^ 


Ax cos Pi -b 


KS 

7% 


cos(pa - 




Vd 

Vu 


1 



(37) 


The scaling effect on the potential corrections can be neglected as being of higher order. Then comes the scaling 
of the mass-matrix: 


As 


^^ Aa cos Pi -b - j = cos(pA - ( fi ^) + 


Vd ^ 

Vu 

1 ^ 


(38) 


so that we can extract the DR squared mass for the physical state via a /3-angle rotation. Finally, the Z-factors 
have to be subtracted from the pole corrections, since they have been accounted for elsewhere. This provides: 


^a\z 


Xs 


\J ZsZh^Zhj^ 


Ax cos Pi 


KS 

7^ 


cos(pa - Pk) + 


p„-i-Pd 

VuVd 


1 -b ^poie -b cos2 /3 (Zjj^ -1) -b sin^ /3 ( Zh^ - 1) 

(39) 

differ by a factor 1 - ^ -b {6Zh^-6ZhJ. 
This explains the mismatch, reaching the order of one-loop effects, that is 0(1%) here. In particular, the steeper 


Expanding the Z-factors as Z. = l-bi5Z., we see that Eqs 


36 


and 


39 
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m|, (GeV) Precision 2 



(GeV) Precision 2 



Figure 2: Details on the Higgs masses of Fig{^ On the upper left-hand quadrant, we show the Higgs masses close 
to ~ 125 GeV for precision 0 (green lines), precision 2 (blue lines) and for our implementations (red dots). The 
plot on the upper right-hand side compares our results (dots) for the ‘heavy’ masses with those of precision 0 
(green lines). The same exercise is carried out in the lower left-hand corner for precision 2 (blue lines). In the 
lower right-hand quadrant, we alter our implementation of the Higgs mass corrections so that all p^-dependent 
terms are taken into account as pole-corrections only (so that the wave-scaling factors are set to 1 ): the results 
are displayed as khaki dots and compared to the masses obtained in the procedure of precision 2 . 


apparent slope with varying A, in Figj^ is largely driven by the Zg factor. In principle, the approach including 
the wave-function scaling factors is the most refined among the two methods, hence should be prefered. On the 
other hand, our choice of setting the Z-factors at a low-value of the external momentum, = 125 GeV, is not 
optimized for heavy states. In any case, a 1% effect should not be taken too seriously in view of the various 
additional sources of uncertainty (parametric errors, running, etc.). 

We then consider a second example with A = 0.7, n = 0.1, tan/3 = 2, ~ 2.33 Ma + 20.45 GeV, Ma € 

[0.3,3] TeV, = -50 GeV, = 0.5 TeV, = 1.5 TeV, = 110 GeV, 2Mi = Ma = 150 GeV, 

M 3 = 1.5 TeV At^b.T = —0.1 TeV. The results are displayed in Figj^ This region of the parameter space highlights 
another effect in the NMSSM Higgs sector, namely the large contribution of F-terms to the mass of the SM-like 
state for large A and low tan/3. Indeed, the low value of tan/3, the low mass of the squarks of third generation 
and the moderate trilinear soft terms would result in a Higgs mass below M^ in the MSSM, making this regime 
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Figure 3: Higgs spectrum for A = 0.7, k = 0.1, tan/3 = 2, Ma — 2.33^eft + 20.45 GeV, Ma G [0.3,3] TeV, 
= -50 GeV, 5 = 0.5 TeV, = 1.5 TeV, = 110 GeV, 2Mi = M 2 = 150 GeV, M 3 = 1.5 TeV 
^t,b,T = —0.1 TeV. (Comparison of our results (dots) with precision 0 (upper left-hand plot), precision 2 (upper 
right-hand plot). Focus on the masses close to 125 GeV (bottom left-hand plot). Finally, results obtained with 
HiggsBounds and HiggsSignals. 


incompatible with LEP limits and the LHC measurement. In the NMSSM however, we observe that the mass 
of the SM-like state remains above 120 GeV: this is a consequence of the specific tree-level contributions to the 
Higgs mass matrices, associated with A. Comparison of our results with the masses obtained with NMSSMTools for 
precision settings 0 and 2 again show that our calculation is typically closer to precision 2 , although the differences 
are larger than in the previous scan (about 1 GeV for the two light CP-even states, as can be observed on the plot 
on the lower left-hand corner). We also display the output of HiggsBounds and HiggsSignals for our results (plot 
on the lower right-hand side): HiggsBounds exclusions apply e.g. in the presence of very light Higgs-states with 
non-vanishing doublet composition. The t^st of HiggsSignals provides values down to ~ 75 - for comparison, 
we obtain ~ 78 in the SM limit - when a light doublet is present close to ^ 125 GeV. 

ii) Higgsino and gaugino masses 

Our implementation of the chargino, neutralino and gluino masses should prove very similar to the original 
subroutines within NMSSMTools in the CP-conserving limit. Nevertheless, small technical differences should be 
noted: 
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Figure 4: Higgsino / gaugino spectrum for A = 0.55, k = 0.45, tan/3 = 12, G [100,1500] GeV, Ma = 1 TeV, 
= -300 GeV, rrif g = 1 TeV, Q = 1.5 TeV, mi ^ = 200 GeV, 2Mi = Ma = M 3/3 = 0.5 TeV, 
= —1.5 TeV. Comparison of our results (dots) with the implementation within NMSSMTools (solid lines; 
note that we actually display the absolute values of the masses). The plot below shows the ratio between our 
results for neutralino masses and those delivered by NMSSMTools. 


• we take into account the Higgs-higgsino-singlino couplings which had been neglected in NMSSMTools: this 
results in additional corrections to the higgsino and singlino masses; 

• similarly, bino and winos are not assumed degenerate in the calculation of loop corrections to the higgsino 
masses; 

• all masses are chosen real and positive: this is possible since the diagonalizing matrices are complex. The 
convention in NMSSMTools consisted in keeping these matrices real, so that some masses could take negative 
values. 

We consider the following region in the NMSSM parameter space: A = 0.55, k = 0.45, tan/3 = 12, G 
[100,1500] GeV, Ma = 1 TeV, = -300 GeV, uif g = 1 TeV, = 1.5 TeV, = 200 GeV, 2Mi = 

M 2 = M 3/3 = 0.5 TeV, At^b.T = —1.5 TeV. The masses of the higgsinos and gauginos are shown in Figj^ The 
scan over drives a significant variation of the higgsino masses, while the gaugino masses remain essentially 
constant. Once again, the masses obtained with the original routine of NMSSMTools are depicted with a solid line, 
whereas our results appear as dots: the general features are identical. More quantitatively, the main deviation 
reaches ^ 3% at the level of the neutralino masses: it originates from the corrections to the singlino mass, which 
were neglected in NMSSMTools. 

For the rest of the spectrum, e.g. the sfermion masses, our calculation reduces, in the GP-conserving limit, 
to the original implementation within NMSSMTools. Therefore, we will not push the comparison in this limit any 
further. 
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4.2 CP-violating case 


CP-violation could induce several phenomenological effects at colliders. The most immediate one would be the 
measurement of EDM’s. The absence of any hint in corresponding searches thus places stringent limits on new- 
physics phases. Note however that, at one loop order, these effects are essentially driven by the gaugino phases. In 
other words, new-physics phases associated to the Higgs sector or the third generation sfermions enter the EDM’s 
at the two-loop level only and are thus more loosely constrained. CP-violation could also intervene in rare flavour 
decays and oscillations, which are consistent so far with the SM-interpretation (where only the CKM phase is 
present): such effects have not been included in our study yet and we will not discuss them here. 

i) CP-violating effects in the NMSSM Higgs spectrum 

CP-violation could enter the Higgs sector at tree-level, via a non-vanishing phase ipx — or at the loop-level, 
e.g. via the phases associated to the sfermions of third generation. As a first consequence, the neutral Higgs states 
would become scalar / pseudoscalar admixtures, which affects their couplings to SM particles: for doublet states, 
the pseudoscalar component does not couple to a ZZ or W~^W~ pair, so that the corresponding decay channels, as 
compared to the fermionic decays, are suppressed / enhanced with respect to the case of pure CP-even / CP-odd 
eigenstates. Other effects can be measured in the fermionic channels, provided, however, that the fermion masses 
are sufficiently large. Therefore the presence of CP-violation in the Higgs sector could be tested in precision 
analyses of the Higgs properties - for the observed or hypothetical new states. Note however that doublet Higgs 
states are typically shielded from CP-violating mixing - consider e.g. the zero-entries in the tree-level mass-matrix 
of Eqj^-, so that only a very high degree of precision in the measurement of the branching ratios would be likely 
to detect the tiny - radiatively-generated - pseudoscalar component of a mostly CP-even state. Moreover, the 
current limits on Higgs searches tend to favour a sizable mass-hierarchy between the SM-like Higgs state and the 
approximately degenerate ‘heavy-doublet’ states. This makes the presence of a pseudoscalar doublet component 
within the observed Higgs state unlikely, as the mixing of this state with the ‘heavy-doublet’ pseudoscalar would 
be suppressed in proportion to the large mass gap. Another test would involve the two ‘heavy-doublet’ neutral 
states, which are generically close in mass, so that their mixing could be significant. Yet, the detection of CP- 
violation there will still require high-precision experiments (and the discovery of these states), due to a typically 
reduced production cross-section - with respect to a SM Higgs boson at the same mass; this is related to the 
mostly iJrf-nature of these states - as well as the opening of many less-controlled decays (e.g. towards new-physics 
states). 

In the NMSSM, another type of CP-violating mixing is allowed: a mostly CP-odd singlet may mix with the 
doublet CP-even states - provided A and k are large and (fx — is non-vanishing - and this effect could be fairly 
important if these states are close in mass. In the following, we focus on the SM-like Higgs state at ^ 125 GeV. Such 
a scenario is studied in Figj^for A = 0.68, k. = 0.1, tan/I = 2, = 635 GeV, Ma = 1.5 TeV, mf ^ - = 0.5 TeV, 

mp E = 1.5 TeV, 2Mi = WI 2 = M 3/3 = 0.5 TeV, At^b,T = —0.1 TeV. CP-violation is induced through variations 
of Pk.: note that this strategy is the safest in view of the EDM’s, as non-vanishing ifx produces direct CP-violation 
in the doublet higgsino sector (as well as in the sfermion sector). In the first column of Figj^ = —100 GeV, and 
the mostly CP-odd state is relatively far in mass (~ 250 GeV): correspondingly, the mixing with the SM-like state 
does not reach 1%. The latter state has a somewhat low mass of ^ 121 GeV which translates into a mediocre fit to 
the LHG-observed signals, hence a high x^-value with HiggsSignals. In the second column, we take A^ = 0 GeV, 
so that the CP-odd singlet is close in mass to the SM-like state: at = 0, the singlet has a mass of about 
~ 115 GeV. Consequently, a significant mixing develops between the two light states as soon as 7 ^ 0, the effect 
reaching the level of 30 to 40%. A consequence is the uplift in mass of the heavier SM-like state so that the 
associated signal gives an improved fit with the LHC data. The column on the right is obtained for A„ = 10 GeV: 
the GP-odd singlet is then somewhat lighter (~ 100 GeV), so that the mixing effect at non-vanishing ip^. remains 
milder than in the previous case, yet generates an uplift of the mass of the SM-like state as well. It is to be 
noted that the mostly CP-odd singlet acquires a CP-even doublet component which reaches 0(10%) (at the level 
of the squared mixing angles): the latter would generate a signal at the O(10%)-level as compared to a SM-like 
state at the same mass - indeed, the production cross-section at colliders is essentially mediated by the doublet 
components. For a state with mass ~ 100 GeV, the corresponding signal could be consistent with the LEP ~ 2.3 cr 
excess in Higgs searches with a bb final state 
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even though the state is dominantly GP-odd. 


Note that the two effects that we highlighted - uplift of the mass of the SM-like state via its mixing with 
the singlet and presence of a ‘miniature’ Higgs boson under 125 GeV - are well-known in the GP-conserving 
NMSSM [^, provided the auxiliary singlet is GP-even. GP-violation extends this possibility to CP-odd singlets. 
Further consequences appear on Figj^at the level of the branching fractions of the Higgs states - we display their 
values for the cc and 77 final states -: similarly to the case where the SM-like Higgs boson mixes with a 
CP-even singlet, the proportions among doublet components and may fluctuate, displacing the branching 
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Figure 5: Characteristics of the light Higgs states for A = 0.68, k = 0.1, tan ,8 = 2, /i,ff = 635 GeV, Ma = 1.5 TeV, 
= -100, 0, 10 GeV, mf ^ - = 0.5 TeV, = 1.5 TeV, 2Mi = Ma = M 3/3 = 0.5 TeV, At,b,r = -0.1 TeV. 

The plots on the first line show how the masses of the Higgs states close to ~ 125 GeV vary with the second 
line displays the magnitude of the scalar components of these states. Characteristics of the lightest state are 
shown in blue and those for the second lightest are shown in green. Finally, the lower series of plots shows how the 
points compare to phenomenological limits: the violet mark indicates that the points are excluded by the test in 
HiggsBounds while the brown mark stands for tensions with the EDM’s. The x^-test of the Higgs data is obtained 
with HiggsSignals and corresponds to the red curve. 

BR(S° —bb) 10 GeV BR(S°^cc) BR( S“ ^ YY) 




Figure 6: Branching ratios of the light Higgs states for = 10 GeV in the scan of Fig[^ 


ratios. However, the main effect in Fig|^ concerns the rates of the lighter singlet state which become dominated by 
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CP-even-like channels - for fermionic final states, rates differ at the radiative level depending on the CP propertj0 
while the fluctuations of the branching fractions of the mostly CP-even doublet are dominated by the variations 
of the associated Higgs mass. 

Disentangling this scenario - where a light mostly CP-odd singlet mixes with the SM-like Higgs boson - from 
the CP-conserving one - where the light singlet-like state is genuinely CP-even - is likely to prove very difficult. 
The reason rests with the observation that the singlets do not lend specific properties to the SM-like Higgs state - 
they simply reduce its total width and might alter its branching ratios at the percent level. Moreover their decays 
are essentially mediated by the doublet component which they acquire in the mixing, i.e. a CP-even one in both 
cases. Typical singlet decays - towards hypothetically lighter singlet states or singlinos - would not necessarily 
help to discriminate among CP-even and CP-odd mixing and would be problematic in terms of compatibility with 
the measured Higgs signals. Indeed, the standard rates would then be suppressed in proportion of the magnitude 
of the unconventional decays. While deviations of the rates of the observed Higgs state from the standard ones 
might be interpreted via such a mixing effect - should such deviations be detected at the LHC or a future linear 
collider -, it is questionable anyway whether the light singlet could be detected - possibly in Higgs-pair production: 
see e.g. in the CP-even case. 

At the outcome of this discussion, we see that, while the CP-violating effects involving singlets in the Higgs 
sector may be larger than in the pure doublet case, they are also more difficult to trace and could be mistaken 
for CP-conserving phenomena. For this reason, it is essential that CP-violation be tested in processes where the 
CP-properties are well-controlled, which brings us back to EDM’s or rare flavour transitions. Spectral effects in 
the Higgs sector are unlikely to allow for discrimination with the CP-conserving case. 


ii) Comparison of the Higgs mass predictions with the existing literature 
We will now compare some of our results with existing analyses in the literature, where CP-violation has been 
considered. Note that, contrarily to the comparison with the calculations in the CP-conserving NMSSMTools, one 
should not expect much more than a qualitative agreement. Indeed, the choice of disparate procedures in different 
tools, e.g. concerning the definition of the input - such as the choice of running Yukawa couplings or that of 
versus cos (p 2 -, are known to lead to sizable deviations, already in the CP-conserving case. The level of 
precision in radiative corrections is also to be considered. 

NMSSMCALC is a public tool computing the Higgs spectrum and decays in the Za-conserving but possibly 
CP-violating NMSSM. The chosen approach is that of a diagrammatic calculation. The level of precision has 


recently been extended to include the dominant two-loop corrections 37 


First, we focus on the results of 37 dealing with CP-violating effects, i.e. essentially Fig.6 and the surrounding 
text in that paper. If we blindly input the parameters given in section 4.1 of this reference into our framewor]|^ 
the spectrum ~ not unexpectedly and already with the CP-conserving NMSSSMTools - turns out to be slightly 
different from the quoted one: in particular, the mostly CP-even and mostly CP-odd singlet states appear with 
masses ~ 108 GeV and ~ 113 GeV respectively. Yet, this discrepancy can be absorbed within a small shift of 
cos(p 2 : using the value 203 GeV, we then recover states at ~ 103 and ~ 128 GeV so that the Higgs spectrum 


then largely coincides with the one provided in 37 


In any case, this manipulation has little effect on the properties of the mostly /i°-state, close to 125 GeV. 
Scanning over the phases ip At , T’/i “ a scan over ip^ in the notations of 37 would correspond to a scan over 


ip\, keeping = ipx, in ours, so that GP-violation does not enter the Higgs sector at tree-level - and 4>m3, we 
obtain the plots of Figj^ On the upper part, we observe that the general dependence of the ‘SM-like’ Higgs mass 
on ipAt and ip^ is largely reminiscent in shape and magnitude of that observed in Fig.6 of 37 . In these two 
cases, CP-violation enters the Higgs sector via radiative corrections, where the leading effect is generated by the 
sfermion corrections. On the other hand, the mass obtained with our code is independent from while such a 
dependence already appears at one-loop in 37 . Note that one does not expect gluino corrections to the Higgs mass 
at one-loop order and it is_thus not surprising that our implementation does not show any variation with ■ The 
corresponding effect in 


37 


is explained there as an artifact of the top-Yukawa DR counterterm-fixing of higher 
order. Note also that the corresponding fluctuations, at the GeV level, are small compared to the uncertainty that 
one naively expects for the Higgs mass (a few GeV). 

In addition, we show the values of the EDM’s that we obtain in these scans. These have been normalized to the 
experimental upper bounds: ^ 1 • 10“^® e cm for the electron 
^ 1.3 • 10“^^ e cm for the Thallium atom 
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38 


- estimate from thorium monoxide experiment 
3.1 • 10“^'^ e cm for the Mercury EDM 40 and ^ 3 • 10“^® 


e cm 


for the neutron 41 . Note that only the central values are displayed, without error bands. The color code is the 
same as in Fig.6 of 37 , i.e. green for the scan on ip^^ red for that on ipAt and blue for the one over (when 


®There is also some difference at tree-level, but the corresponding effect is very small for light fermions. 

®Note that this addresses the Di?-parameters in the reference, since the parameters within NMSSMTools are regarded as DR. 
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Figure 7: Top: phase-dependence of the mostly Higgs state as obtained from section 4.1 of 37 . Below: 


estimates of the EDM’s (normalized to their experimental upper bounds: see text). The color coding follows that 
of Fig.6 of [^: green for a scan over = ipx = red for a scan over pAt and blue for a scan over (j)M 3 - 
When the blue curve does not appear, the reason is that associated values are negligibly small. Note that several 
estimates are employed for the neutron EDM. 


the curve does not appear in the plot, this is because the corresponding values are negligibly small). We see that 
the scan over (p^ may generate tensions with the EDM’s - mostly the electron EDM - when is not trivial. 

We now turn to the one-loop analysis proposed in 42 . We first consider the scenario presented in section 


4.1.1 of this reference, where CP-violation intervenes in the Higgs sector at tree-level via the phase Again, a 
qualitatively close spectrum can be recovered with little alteration of the input proposed in the reference and our 
results are displayed in Fig|^ while small differences appear, both the Higgs masses and the composition of the 
states agree reasonably well with those of 42 . The major source of deviation is associated to the use of different 
input - in NMSSMCALC instead of cos p 2 in our case so that the comparison makes limited sense when p 2 
becomes large (i.e. for ~ t^/8). In the regime considered here, the CP-even and CP-odd singlet states are close 
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Figure 8: Masses and CP-even composition {X([ ^ ^ -|- X^ for the three lightest Higgs states in the 


scenario of section 4.1.1 of 42 

nig (GeV) 


nig (GeV) 



Figure 9: Masses in the scenarios of section 4.1.2 (scan over (p^ = (px = Pf^) and section 4.1.3 (scan over pAt ) 
of [421. 


in mass to the SM-like Higgs boson, so that the non-vanishing generates a substantial mixing of these three 
states. 

[4^ then considers the case where CP-violation is absent in the tree-level Higgs sector, but radiatively generated 
via phases in the supersymmetric spectrum. In the first case (section 4.1.2), the ‘active’ phase is px but the 
condition p^^ = px ensures that no CP-violation enters the Higgs potential at tree-level - we will recycle the 
previous notation p^ for this scenario. In the second case, only the phase pAt is non-trivial. We display our 
results in Figj^and observe that they capture the effects depicted in Fig.5 and Fig.7 of 

Our code is thus able to reproduce the main qualitative features that were observed in the CP-violating case 
by NMSSMCALC analyses. We stress that a more quantitative study would have limited interest, as the divergent 
treatment of the input already generates discrepancies between the CP-conserving NMSSMTools and NMSSMCALC. 
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5 Conclusions 


We have presented a series of Fortran tools extending NMSSMTools to the CP-violating case. Radiative corrections 
to the supersymmetric and Higgs masses are computed at one-loop order. Dominant two-loop effects to the Higgs 
masses are also included in the double-log approximation. Additionally, Higgs couplings and decays, as well as 
top two-body decays and EDM’s are implemented and allow for phenomenological tests of the spectra. We have 
shown that our code compares competitively with existing results, both in the CP-conserving and CP-violating 
cases. The new tools will be made publicly available on the NMSSMTools website 13 in the near future. 

We also highlighted a scenario made possible by CP-violation, where the SM-like Higgs would mix with a 
mostly CP-odd singlet state. The consequences on the Higgs phenomenology are similar to the CP-conserving 
mixing with a light CP-even singlet so that both scenarii should prove difficult to discriminate, unless genuine 
CP-violating effects - e.g. in EDM’s or flavour physics - are discovered simultaneously. 

Einally, we would like to close this discussion with some details concerning the future developments which 
we plan to consider. First, an extension of our tools including Zs-violating terms should raise little difficulty. 
Then, flavour constraints are relevant in the CP-violating NMSSM and we intend to design phenomenological 
tests accordingly. Finally, the dominant two-loop corrections to the Higgs masses will be calculated in a more 
quantitative way. 
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A Reference functions 


To{m^) = w? 


TO 


TT / 2 2 2\ 

> 5 (TOi,TO 2 ,TO 3 ) 

TT / 2 2 2\ 

> 6 (TOi,TO 2 ,TO 3 ) 


1 




m 


m 




1 1 

f m2 


(to2 — M^)^ 


- m 2 In 


m2 


-1 


(m? - ml){ml - ml){ml - ml) 
ml In Too li 

-L m.i 2 


(m2 — m|)(m2 — 7713)2 


m 


In ™_0 

TO ^2 ^ 


-M2 


(m^ — TO2)(m2 — 
In^ -2 


m 


2'i2 


(m2 — M2)3 


1 

(m2 — ml){ml — ml) 


We only consider the finite part of the loop integrals: 


^o(to) = -Jo(to^) 

xm2 + (1 — x){M‘^ — p2) + (1 — a;)2p2 


Bo{p,m,M) = - f In 
Jo 

/•i ; 

Bi{p,m,M)= / (1 —x)ln- 
Jo 




dx 


+ (1 — x){M'^ — p2) + (1 — x)2p2 


■da: 


Finally, we borrow some of the notations of ^f\ : 

Bff{p, to, M) = (p2 _ 77^2 _ _ Ao(to) — Ao{M) 

Bsv{p, TO, M) = —(2p2 + 2m2 — M‘^)Bq{p, m, M) + 24 o(to) — 2Aq{M) (Feynmann gauge) 


B The tree-level masses and couplings 

This appendix provides the reader with a detailed presentation of the tree-level spectrum and couplings of the 
CP-violating, minimal-flavour-violating, R-parity and Z 3 conserving NMSSM. 

B.l Tree-level masses 

Here we derive the tree-level bilinear terms of the lagrangian. For a later application to the Higgs couplings as 
well as to the loop-corrections in the Coleman-Weinberg effective potential, we will try to keep a full dependence 
in the Higgs scalar fields S, = (M+,M°)^ and Hd = {H^,HJ)'^. To evaluate the masses, one of course simply 
needs to replace these fields by their v.e.v.’s. 

B.1.1 SM fermions 

The Higgs-fermion potential reads: 

Vf = - HIul)Y^uI + (FfOrfi - HjuLVcKM)Ydd% + (M^e^ - HjFF)Y,e% + h.c. (40) 
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Focussing on the third generation (and neglecting off-diagonal CKM elements), we may cast under matrix form: 


Vf3 = {uL,u%.,dL,d%)Mq3 + {vL,eL,e'^jf)Mi3 


Mq3 = 


0 

YtK 

0 

-YbHl 


from which we derive the squared-mass matrices: 


Ml3=M\,Mq3 = 

0 

-r/ijOif- - Y.^Hl^H- 


0 -YtH+ 
-YtH- 0 
0 YbH^* 


-Y,H+- 

0 

0 


Mi3 = 


-Y,H^ Y,H^/ 


Y,\m-^ + H+H-) 
0 




Yi\H^,? + Y,^HiH- 


0 -Y^H+ 

0 YrH° 

H°* 0 


-YtY,{H^H+ + H^^H+) 
0 

YtHm\^ + H+H-) _ 


r y:^h+h- -y^h°h+ 

Mf3=Ml,Mi3= -Y^H^^*H- i;2|i?0|2 


F,2(|7j0|2 + i7+i/-)^ 


Replacing the Higgs fields by their v.e.v.’s, one obtains diagonal matrices (-^gS/ and (A^p), with the usual 


relations: ml = Y^v"^, ml = Y^v^, ml = 0, = Y^^vl- 


B.1.2 Electroweak gauge bosons 

From the Higgs kinetic terms, one obtains the Higgs-gauge potential (where we omit the derivative Higgs couplings): 

Vg = 1 ■ 1#^] (|//„|2 + \Ha\‘^) + (44) 


After the helds are rotated to the mass-states, 
^ gB^+g'W^ 


-g'B^ + gW^^ 


; ' ^ ^ ; Z^= ^ 7 ^ ' ; W^ = -^ 

/F+F VF+F 


Wl T iWl 

V2 


we derive: 


Vg = ^{A^.Z^) 


0 0 
0 


'■ ] + \Hl?) + ( ,/(>^>) if-f? I + HjHJ) 

^ \ 32+3,2 2(32+3,2) / 

+ y (l^°P + Wl? + H+H- + H+H-)W;W+^ 


V2^g^+g'^ 


{gA^ - g'Z>^) [W+iH-Hl, - H^HJ) + - H^Hj)] (46) 


This leads to the usual gauge-boson masses: = 0, = \{vl -|- u^), Mf = ^ {vl + v^). 
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B.1.3 Sfermions 

The Higgs-sfermion potential originates from soft, F and D terms: 

Vp =mQQpQL + + mjjDp Dp + m\L\LL + m%Ep Ep 

+ {H°Ul - H+DlV^ppj)YuAuF^^-U^p + (H^Dl - HJULVcKM)YdAdF^^^ 

+ {H^^El - H-Nl)Y,A,E^^^E^p + h.c. 

- HIUl)Yu^ + \Y^VcKMHtU‘p - YdH^^D%\^ + \YuHIUI - YdHJV^^^D^p 
+ \{H^dDL - H^ULVcKM)Yd\' + \iH°EL - HjNL)Ye\' + \YeH^E^p\‘' + \Y,H°E^p\‘' 

+ \XF'^^SH+ - ULVcKMYdD'p - NLYeE'pl'^ + - DLYdD'p - ElY^E'j^^ 


(47) 


+ \XF^^SHI-UlYuU'p? + 


XF^^SH- - DpVlp^Y^U^p 


9' 


9 


HlHu - H\Hd + \qIQl - - l\Ll + 2E^JE 


R 


m-^H^ + Hl-^Hd + QI^Ql + 

O I 

The bilinear sfermion terms can be cast under matrix form: 


VpB{ulU^p,DlD'p) 


Ml 


Ml\u 

Ml 




I FiNlElEl) 


\D^JJ 


■^%N 





(48) 


with the matrix blocks: 

,2 I v2\ ir0|2 


Ml = 


m: 


Q 


^ U P-'ttI 


- g^) (|7J0|2 - |7J0|2) r„ 


Yu - Xe-^‘<^^S*H°*' 


jVcKMYiVlp^H+H- 


+ l(^+5") iH+H- 


^u + YulK 

-HjH-) 


0|2 


0 

Y^H+H- - S^iH+H- - H+H^ 


(49) 


Ml = 


+ Yim + i (^ + 9^) - m) Yd [Ade-^^-^nr - XF^^^SH^u] \ 

Yd [AdF^-^H^, - Xe-^^^S*H°^*] ml + Y^\Hl^ + ^(ji/Op - \Hl^)) 

f V^kmY^VckmH+H- + I (^ - 9^) {H+H- - H+HJ) 0 

I 0 YiH+HJ + [H+H- - H+HJ) 


(50) 


Ml,u = 


-V^kmYuHIHF - Y!vIp^HI*H^ + + H^HJ) -V^^p^Yu [Aue-^^-^R- + 

^ -Yd [AdF^-^HJ + Xe-^^^S*H-] -YdV^c pmYu{HIHF + 

(51) 


Ml = ml - 


- \Hl\l + Y^H+H- + - H+H-) 


Ml = 


_(ml + y, 2 | 7 j 0|2 + =^(|^ 0|2 _ |^ 0 | 2 ) - XF^^^SH^ 

Y, - Xe-^^^S*H^u*] + if - W), 

-^-^{H+HF -H+HJ) 0 \ 

0 Ye^H+HF + ^ {HlHF - H+HF )) 

M%p = (- (if - HTHF + ^^HlHF -Fe [AeF^-^HF + Xe-^^^S*HF]) 


+ 


(52) 

(53) 


(54) 


Moving to the v.e.v.’s, the matrices become block diagonal - each block being associated to a given electric charge 
of the sfermion fields. Under our Minimal Flavour Violation hypothesis the various generations also decouple so 
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that we are left with 2x2 (hermitian) mass-matrices Those can be diagonalized via unitary matrices , 

according to: 


{Ml) = X^^diag{m\,ml^)X^ 


2 1 
= 2 


X^ = 


cosOf — sin0i7’e*‘^^\ 
sin^Fe”®'^^ cosOp ) 


^ (■^f)22 - \/((-^f)ii - (-MlW + 4 (MDf, 


1- (Ml)^^ + \J— (MD^I) + 4 {md^^ 


\ Op = arctan 

= arg[(A' 


[ 2|(Ar|)^j 

*f)i2] 

J 


(55) 


The mass-states are given by Fi = Xf^Fp + Xf^^F^* (where, in our notation 1 o L and 2 o i?). 


B.1.4 Charginos and neutralinos 

The gaugino-higgsino-Higgs potential may also be cast under matrix form: 


Vx = oX" 


M^+~ 
, Xi^o- 


M^-+ M^-o 


M^+o X + ^-c- 


= (-zw , , -iw+, h+, -zb, -zw^, h°, h^, /i°) 


Xiy.O + 


M^o 


M 


- i gHi* 


Xi^o = 


M^-o = 


Xi^+O = 


0 

9 fjO* 

g' zjo * 
~72^d 


gHl* 

0 gH+ 


= Ml+- 


-9_}J+ _ 

F2^d 

-y-.^d 0 

0 

XFv^H+ 

0 

0 gH- 

0 

0 \ 

9_fi- -a-H~ 0 

F2 ^ 72 “ ^ 

0 

XF^^^HJ j 

0 

* 

- 

g' fjO * 

~V2^d 

M2e*'^“2 

* 

1 


9 fjO * 

V2^d 

g ttO * 

72 "“ 

0 

- 

-XF^^S - 


= Ml 


V 0 


9 fjO » 
72^d 


-XF^^S 


0 
0 

,19>\ TiO 

-XF'f’^H?, 


= Mi 


0 -XF^^Hl -XF^^Hl 2 kF'^-S J 


(56) 


i) Charginos 

The 2x2 chargino mass-matrix may be diagonalized via two unitary matrices U and V: 

(M^y-+) = {7^diag(m, ±,m^±)y. To determine rn±, m±, U and V, we consider the hermitian matrices: 

\ X / '' Xl A2 Xl X 2 


(■^x+) = (^x-+y (^x-+) = VUiag{m{^±,m{^±)V 

(X(x-+)^ = C/^cliag(m^±,m2±)C/* 
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which provide: 


2 1 

m‘‘± = 2 

Xi ^ 


2 1 

± = h 

xi 2 



0 v=arctan \ \ x+/.. ^ K\ \ 

(^v = arg[(x2^)j 


9u = arctan 


2\(m^_) I 

\ X / 12 


[(^c/ = -arg[(>[ 2 _^j 

0\ / cos0[/ — sin0c/e*'^^\ A ® ^ cos0y — sin6<ye®‘^''\ 

^ V 0 ly l^sin0[/e-*^^ cos 0,7 ) ’ l^sin^ye"*'^'' cos0y ) 

The choice of phases tpu, ‘fv is a priori arbitrary. We decide to determine them by the requirement that and 
m^±, obtained in the matrix product U* (-^x“+) positive. The associated mass-states are then: 

x’*’ = Vii{-iw~'') + V^ 2 hi = Vinl{-^w^) -I- Viuh'^ ; x“ = U^i{-tw~) + UahJ = Uiw{-iw~) + UidhJ 


ii) Neutralinos 

The 5x5 neutralino mass-matrix is symmetric, hence is diagonalizable via a single unitary matrix N: 
(A^^o) = A^diag(m^o, i = 1,..., 5)N. As before, we first consider the hermitian matrix 

{M\o) = {M^o)^ {^x°) = A^Aiag(m^o,i = 1,..., 5)A 


This hermitian matrix - or equivalently the 10 x 10 symmetric matrix ^ ~ may be diagonalized 

numerically, providing us with m^o, i = 1,..., 5 and a diagonalization matrix Nq. We define N = diag(e*'^’'i, i = 
1,..., 5)Aq, where the phases ip^o , t = 1,..., 5 are determined by the requirement that the masses m^o , i = 1,..., 5 
obtained from the matrix product N* are real and positive. The neutralino mass-states are then defined 

as: 


— Nii(—tb) + Ni2{—tw^) + Ni'^h^ -\- Ni4h^ -\- Ni^h^ = + Niuh^ + Nidh^ + Nigh^ 

B.1.5 Gluinos 

The gluons of course remain massless. Concerning their supersymmetric partners, the gluino bilinear terms read: 

Vg = -M3e^'^^^~ga~ga + h.c. (57) 

so that we define the mass states Ga = —le^’^’^^ga, with mass M 3 . 


B.1.6 Higgs sector 

The tree-level Higgs potential is given in Eq[^ 
i) Minimization Conditions 
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First derivatives of the potential must vanish at the minimum, which provides: 

1 / dVno \ 


= 0 = 


— I'd) 


2 I \2f 2 I 2\ 5 ^+ 5 ^/ 2 2 \ 

+ a (s + u„) -^-(u„ - Urf) 


= 0 = As [Ax sin (/?i + ks sin{ipx - (Pk.)] va 
= 0 = As [Ax sin tpi + ks sin((y9A - Pk)] Vu 


Vu - As [Ax cos Pi + KS C0S((^A - Pk)] Vd 
Vd - As [Ax cos V?1 + KS C0S((/?A - Pk.)] Vu 


2 \ 9/10/72/ 

l/ WgO \ 

2\dhyV2/ 

^ = 0 = [m| + ks{Au. cos P 2 ) + 2ks + A^(u2 + Vd)] s - X [21 acos((?i + 2kscos(v3a - PK)]vuVd 

1 / gVgo \ 

2 \9aS/72/ 

1 / 9Vgo \ 

2\dayV2/ 


1 / 9Vffo \ 


= 0 = A [Aa sin 1^1 — 2ks sin((pA — Pk)] VuVd — ks^A^ sin p 2 


2 \ 900/72/ 

So that one can express certain parameters in terms of the v.e.v.’s: 


rriH^ = As [AaCOSv?! + Kscos{px - Pk)] — - A^(s^ + vj) - ^ t ^ (vl - vj) 

= As [Ax cos Pi + ks cos{px - Pk)] — - A^s^ + v^) + ^ t ^ (vl - vj) 

Vd 4 

rrig = A [^A cos(/?i + 2kscos{px — Pk)] _ ks{Ak cos (^ 2 ) + 2ks — A^(u^ + 

Ax sin pi = -KS sin((^A - Pk) 

Ak sm (^2 = - [^A sm(pi - 2Kssin{px - Pk)] — 5 - = -3A-sm((/?A - Pk) 


a) Charged Higgs 

The 2x2 charged-Higgs bilinear terms read: 




Vh± ■■ 

= {h-,h^ 


9'^ + 9^ 
4 

i]H°]^ 

- wl?) - 

= As [x 4 a cos Pi + ks cos((/?a 

-Vk)] 

^AL-{x^ 

Vu \ 

+ A7|5|2 -s 7 + 

ff'" + 5' 
4 


-]Hl?- 

= ml,,+A2|5p- 

7" + 9 ^ 
4 


- wi?) - 

= As [Aa cos Pi + ks cos{px 

- Vk)] 

Vd \ 


+ g" 
4 

,,'2 I „2 


(58) 




[2H+H- - H+HJ] 
[2H+H- - H+HJ] 




2 0 

2 9^ 

2 


7 + - H+H-] 


r- , 7"+ 7 


[2H+HJ - H+H-] 


+ A 2(|^|2 _ 7) _ 9 ^^ (| h 0|2 _ |^0|2 _ ^2 ^ ^2) ^ ^(|^0|2 _ ^ 2 j 


M 


H±\i2 


= A 


Axe-^^^^S* 


- (A" - + (A" - = Mjj^ I 


= As [Ax cos Pi + KS cos{px - Pk)] - ( A^ - y J VuVd + (A^ - - — X^)H+ 

+ A [^acos(/?i(S'* - s) + kcos(v3a - Pk){S‘^ - 7) - iKsin((^A - Pk){S'^ - sS*)] - (A^ - - VuVd) 


2 9'^+9\„+„- 


V 2 9 \/ jjO TjC 
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Obviously, 


{M‘]j±) = |as [^aCOS (,£)1 + Kscos{ipx - ip^)] “ “ y ) 


1 

Vd 


m^± = 


— sin /3 cos ^0 0 \ /— sin /3 cos /3\ 

cos P sin /3y y J y cos /3 sin /?y 

[Aa costal + Kscos((pA-(/?«)]- I (I’m + ^'d) i tan/? = — (59) 

VuVd V 2 / 1 Vd 


with the Goldstone boson = — sin + cos and the charged Higgs state = cos(3H^ + sm(3H^. 
We will denote the corresponding rotation matrix as follows: 


C — (~ sin P cosP\ 
cos P sin P 


= 


1 ; Hf = 


Hf = H- 


Hi) Neutral Higgs 

The symmetric 6x6 bilinear Higgs matrix Mjjo = I oSi/^V^dl/V 2 ’^ 
following elements: 


includes the 


Mj,o ,,=mj,+X\\S\^ + \H^>\^) + 


0 | 2 ^ , 5 '^ +5^ 


{2Re{H°r + |id°|2 - |if0|2) + - ^-'^-p^H+R- 


^ r A , t , 5 '^ + 5^ 2 , rr+ Tj- 9 '^ - 9^ TJ+ TJ- 

\s [Ha cos Pi + ks cos{px - p^)] — + -^-u„ + -- H+ - - - HJ 


r,l2 




+ - s" + \Hl? - vl) + (2Re(i/°)2 + \Hy - |7J°|2 - + y) 


-^2 — —ARe 




+ 2 I ] Re(H°)Re(id°) 


= -As [Ha cos (,^1 + KS cos{px - Pn)] + 2 ( A^ - ^ 4 ^ ~ y)^®(^« -^d ) 


— A [Ha cos (^iRe(S' — s) + k cos{p}x — :/5K)He(S'^ — s^) + k sin((/9A — V3K)Ini('S')Re(2S' + s)] 

(Re(id°)Re(idO) - v^Vd) - (A^ - ^)Re(iJ+i7-) 


+ 2 (A^-^+^'' 


.21 _ ^2 , , 2 ,, c ,2 , , „o, 2 ^ , 9 '^ ^ {2Re{H')f - lifOp + \HX) + - ^A^h+R- 


22 “ ^Hd + A (|S'| + \Hu\ ) + 


= As [Ha cos Pi + ks cos{px - Pn)] — + -—Hd - - — r-^HpH. 


Vd 


pX^{\Sp - s^ + \Hlp + 


^ a a ^ 

{2Re{Hlf + - \Hy - ivl + vl) 




HO 


13 


= -ARe [(Hac*’’’’! +2Ke*(‘^^-‘^'=)5'*)7d;?J + 2A2Re(S')Re(id°) 

= —Xvd [Ha cos Pi + 2«:s cos((^a — Pk)] + 2A^sv„ 

- A [Ha cospiRe{H)l - Vd) + 2kcos((Pa - PK)RQ{S*Hd - svd) 

+ Ksin((pA — Pk) (lni(id)^)Re(s — 2S) + 2Im(S')Re(7J))))] + 2A^ (Re(S')Re(id°) — sv^ 


Hd^o|23 


= -ARe [(Hab*^! + 2A2Re(5')Re(il)?) 

= —Xvu [Ha cos Pi + 2ks cos{px — Pk)] + 2X^svd 

- X [Ha cos(,?iRe(id° - Vu) + 2kcos{px - (/?„)Re(S'*id° - sv„) 

+ Ksin(v3A - Pk) (lni(i?[()Re(s - 25') + 2Im(5)Re(7d°))] + 2A^ (Re(5)Re(id^) - svd) 
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I 33 =ml + 2KAJie{^^^S) + 2k^ [2Re(5')2 + \Sf] + + \H^f) - 2XKRe 

+ [H+H- + H+HJ] + 2XnRe H+HJ 

= Ks [Af^ cos (p 2 + 4ks] + AAa cos ipi _|_ ^2 \^jj+j^- _|_ JJ +^-j _|_ 2AKRe 

+ 2k cos(^ 2 Re(^ -s) + K (2Re(5)2 + |5|2 _ Ss^)] + [\H°\^ + \HX -vl- 

- 2Xk cos {(fix - Vk)^^{HIH^ - VuVd) - sin {(px - Pk) (Re(i7°)Im(i/^) + Im(7J°)Re(i/^) + 3^^Im(S')^ 
^Uu = ^^^^Re(i/°)Im(7jO) 

A4^o I 24 = -ARe [* [Axe^^^S + 2 )] + 2 (^A^ - Im(i70)Re(iJ°) - (A^ - ^)Im(iJ+i/-) 

= A [(A;^ COS Pi — 2 k cos{px — Pk)^^{S)) Im(S') + k sin((p> — (/ 9 „)Re(S'^ — sS)] 

+ 2 (^A^ - Im(7jO)Re(ilO) - (A^ - ^)lmiH+HJ) 




= -ARe i (^Axe^^^^ + 2Ke*(‘^^-^'‘)S'*) + 2X‘^Re{S)lm{H^) 

= Xnsvd sin(i^A — Pk) + X cos pilm{Hd ) + 2k cos(px — Pk) (Re(5')Im(i/^) — Im(S')Re(7J^)) 

+Ksin((pA - ‘/?K)R-e((2S'* - s)7J^ “ + 2A^Im(_ff°)Re(S') 


M]doU = + l^°n + (2Im(if°)2 + |7J0|2 - |ff0|2) + ^-^h+H- 


_(_ fp- _ ^2 

') + - 3-^H+H2 


/2 2 /2 2 

= As [Aacosv?! +kscos((/?a -<Pk)] — + ^ ^ H + H - - ^ H+HJ 

Vu ^ ^ 

+ + \Hl\^ - vl) + (2Im(i/0)2 + \Hl\^ - \hI\^ -vl + 


A4^o I 45 = -ARe [z S*^)] + 2 (^A^ - Re(i/°)Im(7J°) - (A^ - ^)Im(iJ+i7^-) 

= A [(Aa cos Pi — 2k cos{px — (/?K)Re(S')) Im(S') + k sin((pA — (/ 5 K)Re(S'^ — sS”)] 

+ 2 (^A^ - Re(i?°)Im(ijO) - (A^ - ^)Im(i/+i7,-) 

^ff»l25 = ^^^^Re(i/°)Im(7jO) 

Mldo\^^ = -ARe i (Axe^^^ + 2Ke*(‘^^"‘^»)S'*) 7J° + 2A2Re(S')Im(i7;^) 

= XksVu sin(v3A - Ptt) + A [^a cos (/?iIm(7J°) + 2k cos{px - Pk) (Re(S')Im(7J°) - Im(5')Re(i?°)) 

+/«sin(v3A - t/ 2 K)R-e(( 2 S'* - s)i/° - s?;„)] + 2A^Im(i/°)Re(S') 


A4^o|43 = ARe [Axe^^^^S + Ke^^^-^^^S*^] + 2 (^A^ - Im(i7")Im(i/0) + (A^ - ^-^)Re{H+H-) 

= As [^A cos Pi + KS cos{px — Pn)] 

+ A [Aa cos (^iRe(S' — s) + k cos{px — (/9K)Re(S'^ — s^) + k sin((/?A — </JK)Ini(S')Re(2S' + s)] 
+ 2 (^A2 - lm{Hl)lm{Hl) - (A^ - ^)Re{H+H~) 


34 



+ A2(|5|2 + lijop) + (2Im(i70)2 - \H^J^ + \hX) + - il^H+R- 


/2 2 /2 2 

= As [Aa cos ipi + KS cos((/?A - Vk)] — + ^ ^ ^ H+HJ - ^ ^ H+H~ 

Vd ^ 4 

+ + \H^u\" - vl) + (2Im(i/0)2 + |7J0|2 _ |//0|2 _ 3^2 + ^1) 


■^H»Ii6 


-ARe i + 2A2lm(S')Re(7I°) 

-SXnsvd sin((/9A — V^k) + A [Aa cos (piliii{H^) + 2k cos{(p\ — (p„) (lm(5')Re(i7^) — Re(S')Im(_ff)^)) 

-Ksin((/9A - V 5 K)Re(( 2 S'* + s)H° - Ssu^)] + 2A^Im(S')Re(7J°) 


■^Ho\26 


-ARe i - 2Ke*(^^-‘^'')S'*) i?° + 2\'^lm{S)Re{H^) 

-SXksVu sin((/?A - <Pk) + A [Ax cos ^pilm{H°) + 2k cos((^a - Vk) (lni(S')Re(i7°) - Re(5')Im(i/°)) 

—Ksin((/9A — V5K)Re((25'* -|- s)H^ — 3sz;„)] -I- 2A^Im(5')Re(i7^) 


= 2KAJie [le^^^S] + 4K'^Re{S)lm{S) + 2AKRe - 2A«;Re H+HJ 

= AXKVuVdSm{(px - y^’K.) +2XKsin{ipx- (fift) S^^^VuVd - Re{H°H°) - 2vuVd 

s 

— 2k [A^ cos ip 2 — 2KRe(S')] Im(S') — 2Xk cos(v3a — ^k) [lni(^^°)Re(i?d) + 

- 2 AKRe 

= ARe UAxe“^^ - 2Ke^‘'^^-^’^">S*)Hf\ + 2A2lm(5')Im(i/°) 

= Xvd [^A cos ipi - 2ks cos((px - (/?«)] 

-I- A [Ax cosipiRe{H^ - Vd) - 2kcos((Pa - V 3 /^)Re(S'*i^^ - svd) 

+ Ksin((/?A — ^k) (lm(i/^)Re(2S' -I- s) — 2Im(S')Re(iJ^))] -I- 2A^Im(5')Im(_ff°) 


= -ARe - 2Ke<^^-^'^^S*)H^A + 2XHuv{S)lm{Hf) 

= Xvu [Ax cos ipi - 2 ks cos{ipx - ¥>«)] 

-I- A [Ax cos(piRe{H° - Vu) - 2s;cos((/?a - ip^)Re{S*Hl - svu) 

-{■ Ksin((pA - ‘Pk) (lm(i/°)Re(2S' -I- s) - 2Im(S')Re(7J°))] -|- 2A^Im(5')Im(i7^) 


Mho Igg = ml- 2KA^Re{e^^^S) - 2 k'^ [2Im(5)2 + \S\^] + X‘^i\H^\‘^ + \HX) + 2AKRe 

+ A" [H+H- + H+HJ] - 2XKRe 

= -3KsAf, cos P 2 + X^^A^ [aIa cos pi + 4ks cos{px - Pk,)] + A^ [-ff+ H~ + Hj] - 2AKRe HJ 

- 2kAk, cos(p2Re(5' - s) -I- 2k^ (2Im(S')^ -I- |S'p - s^) -I- A^ -I- -^1.- 'v‘d\ 

+ 2Xk cos {px - pK.)Re{HlH'^^ - VuVd) - sin {px - Pk) (Re{H°)lm{H[l) + lm{H°)Re{H[l) + 3^^Im(S')^ 

As for the case of the charged Higgs, the Goldstone boson G° = — sin /?o° + cos /3a2 can be separated from the 
doublet CP-odd state a° = cos/3a° + sin/3a° by a rotation of angle f3. The remaining (symmetric) 5x5 sub-matrix 

of massive states may be diagonalized (numerically) through an orthogonal matrix : 

{m‘]jo'^ = X^°^diag(TO|o, t = 1,..., 5)X^° (60) 
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The corresponding mass-states are then: 


hi + Xg hi + Xli hi + Xli a° + Xg al ee X^^hl + X^^hl + X^hl + X^y + X^l 


iv) Charged-Neutral Higgs terms 

For completeness we indicate the bilinear terms mixing charged and neutral Higgs states (note that I 




^ -(A" - 

^ -(A^ - 

= HJ + 2X^ReiS)H- 

^ -(A^ - - ^^^HjReiHl) 

^ -(A^ - + ^^^HjReiHl) 

A/("|^+^o = ^ H- +2X'^Re{S)H- 


Iff+aO ~~y ~ ^— 




^ [*A - 2Ke<‘^>--‘^'‘^S*^ HJ + 2XHm{S)H- 


^"\H+aO = ^ -^A" - Y^H-hI - 

^ -*(A2 - + ^-'^^H-lm{Hl) 

^ [-zA if- + 2A2lm(5)ii^ 

B.2 Tree-level Higgs couplings 

Having presented the spectrum and our conventions, we may now turn to the Higgs couplings. 


B.2.1 Higgs-SM fermions 

Employing the Dirac-fermion notation, the Higgs couplings to SM fermions may be cast in the following form 
(with the usual left- and right-handed projectors Pl,r)' 

Vf^![glff'PL + gi”'PR\fS ■ gl^f = [gl*f'fy (61) 

with the (non-vanishing) values of g^^^ : 


S^uu Ti r^R , 1 

9l ~ y 

S°dd Xd T^R , vl ^ 

9l = ^ V^d + yd\ 

S^jee Xe r^R , i 

9l ~ y ^^id] 


gH 


9l ''“ = -nsin/3 


gH--- = -r,sin/3 
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B.2.2 Higgs-gauge 

The situation is unchanged with respect to the CP-conserving case: 

Vg 9 {Sd^^S' - S'8^,3) 


(62) 


r, S'tZZ 

^g^u — gfii 


25 -. 


g'^ + g^ 
' 72 

7 

'72' 

gg' 


S°W+W- it , , vH\ 


HtH 7 ^ _ - ^ 


77^~+7^ 


2ff 


hTh~z 




H^H -77 

25m!^ " = 


9 '^-9^ 

277 ^+ 7 ' 
277 " . 


25 , 


HtH~ZZ (7 -7) 


7^ + 7 

/2 „2'i2 


/2I/ 


— ffAti- 


2(7^+ 7) 


Hf H-W+W~ _ g 

gtiiy = g^vT^Oij 


{vdXl^ + VuX{^) 25772 = J 


77 ^ + 5 ^ 


r " 17 -^ Y"-^ Y"-^ _i_ "Y-^ 

l^id^jd ~ ^iu^ju ~ ^id^jd ^iu^ju. 


„s°H+w- _ 9 

y 


^ [(X« - - (X« + *x7)xfj 


25 




5 ^ 5 ' 


2 gHtZW 


72(7+7": 

99 '^ 


72(7 + 7") 

H+H-Z'f 57(7^ -7) 




7"+ 7 


[u„x£ - UdX)))] 

7 . 


4 7 2-^ _ 9 + ff / yR _|_ vR vR _|_ vI vI I '\ 

9l^y r) -^iu-^ju ' -^id-^jd ' -^iu^ju) 


S°S'^W+W ^ v« I Y« vfl I I '1 

QfJ^y 2 “I” ^iu^ju ' ^id-^jd ' ^iu^ju) 

s«h;w-^ 


gA»^g^7 

72(7"+g": 


-^iu-^ju ' -^id-^ jd ' ^ ju J 

zR I ^ \ vO { vR ^ vl \ x^Cl 


2g, 


s°H+w-z _ 9tii'9'^9 

flU 


72 ( 7 " + 7) 


[(x« + *x7)xf„ - (x« - 
[(x« + ixi)xg - (x« - 


B. 2.3 Higgs-sfermions 

The Higgs-sfermion vertices read: 


Vf 9 SF*F' + SS'F*F' 


(63) 


with: 


gSOutu, =72 


+ 1 - 5 ^ ) (^;„X" - UrfX") 


{xH) 


+ '\/2 




/2 


- ^(u„x« - vdX^^) 


x^R (xYr) 


+ ^ {Xg + iXi) - (s(X« - zXf^) + - *X,4))] 




+ ^ {Xg - zXi) - (s(X« + zXf^) + + *X,4))] (xfj,) 


gS^Dto, ^72 


o'2 


(^/l) 






+ 

+ ^ (X« + zXl^) - (s(X« - zXi) + u„(X« - *X,7)] 

+ ^ (X« - zXl,) - (s(X« + zXl) + v^{X^, + zX^))] X^^ {X^zY 


37 



gS^EtE, ^^2 




+ 




-g' ^ + g" 

4 

J‘2 


(vuX^^ VdXj^^) 


X^L {X, 


YgVdX^ + ^(vuX^ — VdX^) X^j^(X^j^) 


-X [A^e^^PA. {xSd + iXl^) - (s(X« - + z;„(X« - (X^f) 




^ (X^d - ^X^d) - (s(X« + + z;„(X« + (X,^) 


^ ( 9 “ + 9 “) [.uX9 - V,XS\ 


gHtU-D, 


ud 

CKM 


Yu -y) +[yI-\] VdXg Xg^ {XHY + Y^Yd KxS + VdXgY Xgj, (X^nY 


+y„ [A^e^‘^-^Xg + Ase-^‘^^Xg] Xg^ (Xg^* + Yd [Ade-^‘^-^Xg + Ase^‘^^Xg] Xg^ (X,^)*} 

= (g^~^‘^^y 




y2_1_ 


^dXj^d 2 ^uXi^ 


{XYY* + Ye [AeC-^^-^Xg, + Ase*^^x£] (X,^)‘ 




9 


S°S°UlUi 


.'2 


Yg (X«X« + Xix/J + - g2 j (^xY^xY + Xix/, - X«X« - X,^X/,) 

n 


xY, (xYY 


-\x2 ( _i_ \ ^ _i_ ^ 

^u K^iu^ju^ ^iu^ju) '^y^iu^ju^ ^iu^ju~ ^id^jd~ ^id^jd) ^kR\^lR 


1 

+ 2 

+ [{XY - ixDixg - ixY) + (Xg - zXjJixY - ixg)] xYn (xYlY 


4 

Y,X 


+ [{xY + txYYxgd + zxjd) + (Xg + ixgYxYd + txg)] xYl (xYr) 


9 


S?S°DlDi _ 


1 


Yl {xYdXY, + xgxg) + U3-y+gA (X«x« + xgxY - xYdXf, - XgX^,) 


4 V 3 

9Y^ 
6 


XYl {X, 


IL 


\a2 ( vR vR I \ _Y ^ _i_ ^ 

^d y^id^jd ^id^jd) ~^y^iu^ju ^iu^ju ~ ^id^jd ~ ^id^jd) 


xYr {X, 


IR 


+ [{XY - ^xY)ixY^ - ixY) + (Xg - zxjjixY - zxY)] xYr (xYY 

+ [{XY + zxY){xY^ + zxY) + (Xg + zxgxxY + zxY)] xYr {xYY* 


9 


s°s°e:^Ei 


\x2 ( vR yR I yl yl A 1 - 9 '^ + g" / vR vR I vl \rl vR V'i? vl vl \ I \rE { vE 

^e y^id^jd ^id^jd) ;j y^iu^ju ^iu^ju ~ ^id^jd ~ ^id^jd) \ ^kL 


1 

2 
KA 


4 

,,'2 


ILJ 


\r2 ( yR yR 1 yl yl \ 1 ^ fyRyR\y^y^ y^y^ y^ y^ \ 

^e y^id^jd "T ^id^jd) ~:^y^iu^ju ^iu^ju ~ ^id^jd ~ ^id^jd) 


xYr (X, 


ml 


+ XXe-^^Px _ ^xY){XY. - iXY) + (Xg - zXgXxY - zX/j] X^r (X,i) 

+ [{xY + zxYxxgu + ^xjj + (Xg + zxgxxY + zxY)] xYr (xYY* 


= -A^g'^ + g2) [xYxf^-xYxYY 
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9 


S°H+U^Di _ 


j/ud 

^CKM 

V2 


(Yu - - iXl) + {Yi - + tXi^) 


X^L (X, 


+ Y^Y, [Xfjxf, - zXf,) + Xf^iXg + zX^J] X^j, {xEiY 

+Y^\e-^^^xY{XYs - ^XDxYr {xPrY + YAe^^^XgixYs + tXi)X^R (X^)*} 


= (^gH-s»Dru^,y 


9' 


S^H^NlEi _ 


1 

71 


(n" - Y^x^ixY, + zXl,) - - zXl) 


{X. 


IL) 


+Y,\e^^-xY{xY,+zXl) (xYn)*} = (ff 




9 


HtH~UtUi _ 


\/-2\\jrud \2 x^C vO I 

^d \ ^CKM\ ^id^jd 


+ g^]{xY,Xg-XY,xY) 


X^L (X, 


-\^2 \^C Y'C' ^ ( vC vC vC yC \ 

^u^iu^ju ^K^iu^ju ~ ^id^jd) 


X^R (X, 


IRJ 


9 




.'2 


y2\-\/ud \2 yU yU i f ^_ '2 \ ( yC yC V^' ^ 

^u\^CKM\ -^iu^ju' A \ O ^ I \^iu-^ju -^id-^jd) 


rC yC's 


XYl {Xi 


-y2 yC yC i ^_ ( y(J y(J yU y(J \ 

^d^id^jd ~^\^iu^ju ~ ^id^jd) 


C yC 


^CyC^ 


x^R (xPr) 


H+H-ElEi _ 9 '^ + 9^ [vCyC _ yCyC] yB ( vEY 

9 ^ Y^iu^ju ^id^jdj -^kL \^IL) 


\r2 yC yC I 9 _/ x^C yC y(J y(J \ 

^id^jd ~7r\^iu^ju ~ ^id^jd) 


x^R (xYr) 


H. _ \r2 yC yC j_ ( _ /2 , 2'\ [ yC yC _ yC yC 

9 -^id-^jd A \ 9 9 } \_-^iu-^ju -^id^jd. 


B. 2.4 Higgs-charginos / neutralinos 

As for the Higgs-fermion couplings: 


s°x+kxt 


1 


9l 

9l 


[gixYu - ^xLmy;, + - ixl^myCi + + ixDut^vcY 


S°x\x1 _ 1 f ff' 


■^2 I ^ ^ ^ibXD] 

- ^ [(x« - ixk) {n*^n;^ + - (x« - ixl,) {ni^n:, + y^xi,)] 

+ [(X« + zXl) {Nt,Nt^ + N^Md) + {X^a + iXl^) {N^Nl + NtMu) 

+ (X« + iXl) + NlNl^)] 


+ + iXl) [Nt,Nl + N^N, 


ks\ 


^ _ y_xCN*^u*^ + gXYa [nI^U:, - + Ae*^^x£iV*,C/* 

gHrx\xt =y^xY,Nl,V:,+gXY^ (^N:y* + + X^^-xYMsViI 
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B.2.5 Higgs-to-Higgs couplings 

From the tree-level potential of Eq|10[ one may derive the trilinear and quartic Higgs couplings: 


go so so 1 (g'^ + g\ 

/tjSuuu I -rjAuuu YjS udd I fYjSddd , -rjAddd -rjSduu 

9 +n„. -n.,, -n.„ ) + v,{n.„ +n.„ -n,„ -n,„ )] 

I \ 2 r „ /TT'S^ suu I ttA suu i ttS sdd i yrA sdd\ i /TT’S' udd , VfA udd i ttS uss i ttA uss\ 

-h A [5 [llijk -T i-i-ijk "T i-i-ijk "T ^Hjk ) -rVu ^Hjk ^hjk ^Hjk j 


I \ 2 r „ /TT'S^ suu I ttA suit i ttS sdd i ttA sdd\ i /TT’S' udd , VfA udd i ttS tiss i ttA uss\ 

-h A [5 [llijk -T i-i-ijk "T i-i-ijk "T ^Hjk ) -rVu {^Hjk ^Hjk ^hjk ^Hjk j 

+Vd (nf^.f “ + n^.f “ + f * + n^-f ^)] - aH;, cos<^1 ^ 

- a«cos(:^a - cp,) [ 2 s (nf^.f ^ + v, (nf^r* - + 2 n^.r* 

+vu (ng.f^ + 2n^^-f^)] + cos (^2 [nf,r^ - sn^-r*] + 2As [nf^- + 

+ A«sin(^;, - [s + 3n[/“‘^) + {u^r - 2n^r'* - n(,r' 

I fYjPdss qttP sds tt/ dss'i 1 ^u'^d fOTjPsSS Tjl SSS 

-\-Vu y^Hjk ~ ^^Hjk ~ ^Hjk ) ~ ^Hjk 


where: 


N / i,j,k 

(no:-: 


vR vR vR I vR vR vR < vR vR vR < vR vR vR 
-^ia^jb^kc + -^ib^jc^ka + ^ic^ja^kb + ^ia^jc^kb 

+ ^Jl 


^R \rR I vR vR I \rR \rR R 

^jc^kb + ^ic^jb^ka + ^ib-^ja^kc 

{XiXl + XlXib) + 

— (yR yR 1 ~v^P vR\ 1 

V^j6^/cc ^jc^kb) ^ja 

=x! 


yR vR _i_ vR vR \ _i_ f vP vP 

^jb^kc ^jc^kb) ^ja 

+ xl,xlxi^ + xlx^xl, + 


I y/ 

'jc ' -^ic^ 1 

_L vPvP\ _L 
•fcc ^fec 

r _r_ T __ 


_L vP\ 

^ka K^ib^jc ^ic^jb) 

Y^I _i_ ‘v^/ Y'/ _i_ vI vl vI 

ia^'^jc^kb ^ic^jb^ka ^ib^ja^kc 




2„vR , 3'^+5^ 


2A^sX« + 
2A2 sX« + 


_ n'2 

01 _L ^ ^ yR 
-|- VdX^^ 


2 _ /2 /2 I 2 

2„^fl , 5 5 + 


VdXl 


vC vO 
-^ju-^ku 


vC vC 
^jd^kd 


A {Ax costpi -f 2Kscos{ipx - <Pk)) X:[l + 3\Kssm{(px - ^k)X-^ - ^A^ - (vuX^ + VdX^) 

Xks sin((^A - Vk)X^s + ^ (^A cos ipi - 2ks cos((^a - “ y ^ {'^u^ld + VdX-^) 

X I 


g 


30303 ^ 3 ? 


g'^ + 5 " 


[ttS uuuu I TTii UddU OTT*^ uudd , tt 
l^Hjkl ~ ^^Hikl “T 


tS dddd 
^ijkl 


t5 uudd 
^ijkl 


P uuuu _|_ Y^P dddd _ 


ijkl 


^ijkl 


uu P uu _j_ (2y^ 3 dd P dd _ ^'*^^1 

xjkl ijkL xjkl xjkl J 


-f H: 


ijkl 


+ n: 


ijkl 


-hH 


ijkl 


^ + H 


ijkl 


-HH 


P ddss 
ijkl 


+ A^ Aijkl 

- 2Akcos(95a - g^K) 

+ 2AKsin((^A - <^«) - 2ng.;ff«^^ + 


I ttS izu P dd I ttS dd P uu , ttS uu P ss i ttS dd P ss i ttS ss P uu , ttS ss P ddl 

'P^Hjkl ^Hjkl ^Hjkl ^Hjkl ^hjkl J 

I ttP ssud ttS ss P ud ttS ud P ss , OTT'S' su P sd 1 OTT'S' sd P swi 
^Hjkl ~ ^hjkl ~ ^Hjkl ^^Hjkl ^^Hjkl \ 


I 2 FtT'S' ssss I ttP ssss i oTT'S' -P -ssI I 

p [^hjki 'P ^hjki 'P ^^hjki \ f 
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where: 


ttS abed 
^Hjkl — 

\ '' 

/ j '^(T{i)a^(T{j)b^cr{k)c^(T{l)d ’ 

YfP Clbcd _ \ '' \rl -yl yl yl 

^^ijkl / j '^(T{i)a-^a‘{j)b^<T(k)c^(T(l)d 


( 7 ^ iS'4 

crGS 4 

ttS ab P cd 
^Hjkl 

_ \ '' 

/ j '^(T{i)a'^(T{j)b-^(T{k)c-^(T{l)d 

. Tj8aPbcd \ '' yR yl yl yl 

’ ^^^3kl — 2-^ ^(j{i)a^a{3)h^<y{k)c^<yg 


c^e 1S4 

<7^84 

ttS abc P d 
^Hjkl 

_ \ '' 

/ j '^(T{i)a'^(T{j)b-^(T{k)c-^(T{l)d 



<7^.84 



o'2 


H, _ y _ r vR vR I 1 f V’C' vC vC vC\ 

9 ^ l^iu^ju'^iu^ju -^id^jd ^id^jdj \^ku^lu ^kd^ld) 

I 9 f vR vR _i_ _i_ _i_ 1 f vC* _i_ vC vC\ 

■r l^iu^ju “T ^iu^ju ^id^jd ^id^jd\ \^ku^lu ^kd^ld) 

+ {X^uXg + xgxg) 


-jIA^- 


r _i_ 1 f Y'C’ _i_ Y'^^ 

^id^ju ~ ^iu^jd ~ ^id^juj [^ku^ld ^kd^lu) 


I ^ / \ 2 y \ r _i_ _i_ _i_ 1 f Y'C’ vC' wC* Y'C‘'\ 

2 \ ^ / ^id^ju "T ^iu^jd + ~ ^kd^lu) 

+ A«: [cos((pA - + sin(^A - V.) {X^s^l + X^X^,)] {X^^Xg + Xg^Xg) 

+ iXk [sin((^A - <^«) {Xgxg, - Xgxg) - cos(^A - ^«) {xgxg + Xgxg,)] {xg^Xg - Xg^Xgg) 


B.3 Other couplings 

B.3.1 Chargino - Sfermion - SM fermion 

For each fermion / sfermion generation (with the convention of Eq|61[) : 


9l 


= Y^Xgj^-gXg^Vl 


u 


JL' iw 


^ T/' vD tt*- vD rr* 

9l XdXj j^U 9Xj]^Ui^ 

7V*x+e _ Tr* 

9l 9^iw 

Xi >A g ^ ^ 

9l ^e^jR^id 9^jL^iw 


UjXt 

9r 

^jXi 

9r 


= YdXf^Uu 
' = Y^Xg^V,, 


7V*x+e -rr 

9r — ^eUic 

%* x,“ V 
9r = 0 


B.3.2 Neutralino - Sfermion - SM fermion 

For each fermion / sfermion generation: 


9l 

b 

9l 


xN ^ -Y^Xg^Ng - ^ + gX* 


= -YdXg^Xl 


-±f^Xg-gX, 


72 V3 




9l 


7*X?e ^ _y 1 




jL 


vD 
iw I '^jL 


jRXg + ^ { 9 'N*, + gX:j Xf^ 


u^x.u ^ ^ 272^1 

„D*X°d _ ^b XT JvD XT 

9r — ^d^jL^id ^ 9 ^jR^ib 
9r^''' = 0 

= -YeXgXu - V2g'xfjiX,t, 


7*x7 

9r 


B.3.3 Chargino and Neutralino gauge couplings 

Using the notation: 

V, 3 Vjr \gl"‘PL+9l."‘Pii] /' 


9r"=(9F')' 




9r — 


(64) 
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the chargino and neutralino gauge couplings may be written: 


7X, xt 


9l 


99' 


IXi xt 


\/p~+^ 


9r 


99' 


\/g'^ + 


91^'"'^ = ^9^+9'^ {v^n,v;^ + \v^uV*^ - 


9r ' ^ — V9"^ + 9'^ { UiwU*^ + -UidU*^ — s'^Si 


^ix, ^ yy 


x?x“ _ +; 


9r 


\/9^ + 


9‘ 


{N,^N:, - N^dNtd) 


w-xUt ^ g ( ^N^uV*^ - N^^V* 


w-x°xi 


1 

V2^ 


9r - =-9{ -y^UuNtd + U,^N* 


C Radiative corrections to the supersymmetric spectrum 

C.l Electroweak gauginos and higgsinos 


We follow the approach of 22 and consider the loops involving sfermions / fermions, higgsinos or gauginos / 
Higgs and gauge bosons in the self energies of the gauginos and higgsinos, under the assumption that the gauge 
eigenstates are approximately mass states. Taking the complex phases cjiMi and v^a.k into account, we find the 
following corrections to the gaugino and higgsino masses: 

{1 IHi (Ml, 0, MqJ + 9Hi (Ml, 0, Mi J + Bi (Ml, Ai, M^) + Hi (Mi ,g,Mz) 

^ sin2/3cos((()Mi + ‘fx) [Ho(Mi,/r, M^) - Ho(Mi,/i, M^)] | 


/AMi\ 

1 ^ 

V Ml ) 

167r2 1 

fAM2\ 

, 9^ /, 

\ M 2 ) 

' 167r2 1' 


+ 


+ ^ sin2/3cos((/>M2 + [Ho(M 2 ,/r, M^) - Ho(M 2 ,/r, M^)] 

M2 




(t)- * 


4i3i(M2, M 2 , Mw) — 8i?o(M2, M 2 , Mw) | 
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We took over the notations of 22 to designate the approximate masses of the particles in the loops; note that 
pi = Xs stands for the doublet higgsino mass, mg = 2ks for the singlino one and m/^o ,30 for the singlet (pseudo)scalar 
masses. 
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C.2 Sfermions 


0{as) corrections to the squark masses are generated by gluon / squark and gluino / quark loops. Another source 
are squark self-couplings, as these receive a contribution from the SU{3)c D-term. We use the following expressions 
to correct the squark squared masses 
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rn^QBi{mQ,mg,mg) - Aoirug) - ml il + 2^Re{e j Bo(wq, m,) 


We recover the results of 43 in the CP-conserving limit. 

C.3 Gluino 


We follow 22 to include the 0{as) corrections to the gluino mass: these involve the gluon /gluino and the quark 
/ squark loops. The latter depend on squark and gluino phases via the quark / squark / gluino couplings. We 
obtain: 
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Note that the one-loop corrections to the NMSSM spectrum have also been presented in 44 


D Radiative corrections to the Higgs spectrum 

D.l Wave-function renormalization 

We summarize the discussion of section [3. 3. 1[ Remember that = 125 GeV replaces the external momentum. 
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(65) 


Zs = 1 + ^ + K‘^Bo{fiH,ms,mg)} 

D.2 One-loop contributions to the effective potential 

D.2.1 SM-fermions 

The squared-bilinear matrices of the fermions of third generation have been provided in Eq |42[ |43[ One observes 
that they split into (at-most) 2x2 blocks corresponding to the left-handed quark fields, the right-handed ones, 
the left-handed leptons and the right-handed one. The following eigenvalues can be derived: 
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One can then deduce the following (non-vanishing) contributions to the Higgs mass matrices (with = 3 the 
number of colors): 
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We also note the following contributions to the trilinear Higgs couplings: 
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D.2.2 Electroweak gauge bosons 

Using Eq |46[ we obtain the following contribution to the effective potential (where we have performed an expansion 
in terms of the charged-Higgs fields): 
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One derives the following (non-vanishing) contributions to the mass-matrices: 
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Similarly, one can derive the corrections to the trilinear Higgs couplings: 
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D.2.3 Sfermions 


Considering the sfermion mass-matrices of Eq 49 50 51 52 53 54 and focussing on the neutral-Higgs dependence 
first, one obtains decoupling 2x2 blocks - 1 x 1 in the case of the sneutrinos - so that eigenvalues may be 

, m = 1,2. Corresponding contributions to the neutral Higgs 


expressed as m|, = | Tp{H°) + (—l)"*^i?|,(i 70 ) 


mass-matrix thus read (we denote as Eijk the coefficients coming from the tadpole equations - see Eq 32): 
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Additionally, the sneutrino mass reads: ^ (1-^°^ ~ l^dH- 

To derive the corrections to the charged-Higgs masses and Higgs-to-Higgs couplings, one is confronted to the 
task of diagonalizing the matrix-system of Eq|^ This can be performed perturbatively, as an 

expansion in the Higgs-doublet fields, which amounts to a series in We confine to a precision of order O (v^) 

at the level of the masses, which means that we compute the potential up to terms of i7^-order {H standing for 
any Higgs-doublet field) and freeze singlet fields to their v.e.v. s for terms of H'^-ordei (they are kept explicitly for 
terms of lower order in the expansion). The ensu ing corrections to the Higgs potential can be matched onto Eq 33 
and we may then use the results of section E.3 e.g. for the charged-Higgs mass (Eq 67) or the Higgs couplings. 
For squarks of each generation (note that we neglect the Yukawa couplings of the two nrst families): 
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D.2.4 Charginos and neutralinos 


To include the chargino and neutralino contributions to the effective Higgs potential, we turn exclusively to the 
method that we have just presented in the case of the sfermions. In other words, we diagonalize the matrix system 
of Eq 56 perturbatively, in an expansion of doublet Higgs-fields and we match the ensuing potential to the form 
One can then work out the contributions to the Higgs mass-matrix and couplings. Note, however, that 


33 


of Eq 

instead of diagonalizing directly the 9x9 (squared) bilinear matrix of Eq 56 it is easier to consider the dependence 
on neutral Higgs fields only (that is replacing charged fields by 0), as the corresponding matrix then splits into 
various blocks. All the couplings of Eq|^ can be identified from the neutral potential, with the exception of 
A 3 _ 4 , which only appear in terms of the sum A 3 -I- A 4 . It is a straightforward task, however, to compute A 4 in a 
second step, from the charged couplings of the full potential. Another remark accompanies the observation that, as 
higgsino and singlino masses depend on the singlet Higgs field, the coefficients in our matching procedure depend 
on the singlet fields as well (those would correspond to operators of dimension > 4), which leads to additional (but 
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straightforward) terms, with respect to the results of section E.3 This 5'-dependence can be neglected for terms 
of order H^, as keeping it would produce terms of higher order in -rrl—. For and M?, the ^-dependence is 

“sUSY 

largely absorbed (at least at leading logarithmic order) by the coefficients Ap and Ap. We obtain: 
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-45^A^ [(M| + 4KVn A^|5|^ +4KA^M2Re J•6(A^|5|^ M2^ 4K^|5'|^)j 


5 A 5 = 


-2 

32712 


+ 16K^A®S'*'‘e^“^'^^■'^"'J■7(4K^|S|^A^|S|^) + 2g'^/MlM2A^S^e‘^'^“l+'^"2+^'^^^J•6(A^|5'|^Ml^M|)| 


5Xe 7 = tjtA |-!7"‘MlASe‘^'^“l+‘^^'J■3(A^|S|^Mf) + /)MiAS'e‘'‘^"i+'^^’J'3(Mf, A^l"^) 

oiTT^ L 

- 3/M2A5'e'''^"2+'^A)_7r3(;^2|_g,|2^^2^ _ ^2j.^/2 3^2^^^^g,g.{VM2+7>A)_7r3(^2^;^2|_g,|2^ 

+ 8KA®S*^e*'‘^^“‘^'‘' [J•3(AVl^4KVl'^) +•7^3(4/^^|S'|^A^|S|^)] 

-g'^g^XS [Mie4'^“i+'^^) + M2e'''^"2+^^^] J•5(A^|S|^ Mi, M|) 

- g'* {Ml + A^n MiASe*^'^“i+‘^^> J'7(Mf, A^ISI^) - Zg'^ (M| + A^H TjiUl, A^ISI^) 

+ 8«A® (aV4k'') \S\^S*''e^'^^-'<='^'^T7{4.n\S\^,X^\S\^) 

-g^g^XS ]^{mI + X‘^\S\'^) Mae'^^^^+v^A) (m| + A^IS]'') Mi TeiX'^lS]'', Ml, Mi)'j 


D.2.5 Higgs-to-Higgs contributions 

Instead of diagonalizing the Higgs bilinear terms, we compute the Higgs self-energy (H) and tadpole (T) diagrams 
mediated by gauge and Higgs particles in the Feynmann gauge. We then set the external momentum to zero to 
determine the potential contributions to the Higgs mass matrices and subtract the pure gauge effects in the Landau 
gauge - from the results in appendix |D.2.2[ Finally we identify the corrections to the Za-conserving parameters 
of Eqj^ 

i) Pure-gauge contributions to the Higgs self-energies and tadpoles 
For external neutral Higgs: 


n^oso(p") = Ti^ li \g^Ml^Bo{p,Mw,Mw) + ^-^MlBo{p,Mz.Mz) 

^ j ^ L ^ 


-b2 


rpV _ 1 3 

nO — 1 c—O n 


's? 


16712 2 


sin^ P cos 2 pX^^Xf^ + sin p cos -b ) 

g^A,{Mw) + ^^A,{Mz)] + X^Xj^ + X^^X^^ } 


g^Ao{Mw) + ^^Ao{Mz) [vuXg + v^Xl 


For an external charged Higgs (we consider only the physical state): 

2 


H 


V 

H+H 


-{P^) = 


IGtt^ 


g^Ao{Mw) + ' 9 V ,2 \^ ^o{Mz) 

^{g +g) 
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ii) Higgs / gauge diagrams 

The neutral self energy receives contributions from hybrid Higgs (Goldstone) / vector diagrams: 


n 


S?S9 


(/) = 


327r2 


/2 _|_ 2 

g^Bsv{p, Mw, Mw) + - — 7 r^—Bsv{p, Mz, Mz) 


+ 


X [sin" pxtxfi, + cos" IdX^^Xf^ + sin/3 cosp[Xgxf^ + )] 

g^Bsv{p, Mw, Mw) [sin" jSXl^Xj,, + cos" pXiaXjd. - sin/3cosd(-^/u-’f/d + 

+ 5 "Bsv(p, mj^±, Mw) [cos" p{xgxf,, + + sin" PiX^^Xf^ + 

+ sin p cos PiXtxf^ + x 5 x/„ - Xf„X/, - Xf,X/j] 


+ ^^ Bsv(p,mgo,Mz)y. 

k=l 


( ■yr/ I V-^ \ / V"-^ v/’-R V-^ I V"^ V-^ vR \ 

^iu^ku ~ ^iu^ku "r ^id^kd ~ ^id^kd j \^ju^ku ~ ^ju^ku ~x y^jd^kd ~ ^jd^kd j 

For the charged Higgs self-energy: 


^2 5 


j -g^r^-^^sv{.P,mjj±,Q) + 2^^,2 Bsv{p, m^i ,Mz) + ^^Bsvip^mgo^, Mw) 

X [cos" /3(Xf„" + Xil) + sin" /?(X«/ + xL") - sin/?cos/3(X« X« - X^XL)] 


Hi) Pure Higgs loops 

The loops including only Higgs bosons (including the Goldstone bosons, with mass Mw and Mz) read: 


n|;s»(p") -I5„i 


^ra,n—l m—1 

-1-^ s“s;(s;(. o ^ ^ \ SP A Irr, 'll 

+ - 2_^ g ^ rr, r,g 3 ,, "* i^Q (P, , mgO ) + ^ ^ > 

m,n=l m—1 ) 


^H+H-ip") { Y. 


Hi 


^ m—1 


'^o{mjj±) + ^Ys 

m—1 


s'i 


Ao{mso) 


n=1.2 

■El 

771=1,6 




5 *™ " Ho(p,mso ,TOj:^±) 


rpS _ 

^5° - 


ley^Tr" 




m=l 


771 = 1 


where the Higgs-to-Higgs couplings can be found in appendix |B.2.5[ 
iv) Contributions to the Higgs mass-matrices 

Since we are interested in the contributions from the effective potential, we take the limit p" = 0, which simply 
induces the replacement Bq{p, m, M) — )■ M"). The contribution to the mass matrices of the Higgs states 

reads: 


<5 .. = - [n^+J'"+yO) - Tjo+^ 


5mfj± — 


n^+|^+''( 0 ) - cos" - sin^PTl 


where the coefficients Eijt: are the same as in appendix D.2.3 that is, they correspond to the tadpole coefficients 
of Eq[^ 


v) Reconstruction of the Higgs contributions to the potential 
After subtracting the pure gauge contributions from appendix |D.2.2[ one can reconstruct the Z 3 -conserving pa¬ 
rameters of the potential of Eq|^ induced by Higgs corrections. We start by rotating away the neutral Goldstone 
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boson, obtaining thus a 5 x 5 matrix for the Higgs corrections to the neutral Higgs mass-matrix: 6 (M^q) . We 

\ / ij 


then employ the method that was outlined in 28 


SAscos,p, = -i + ^="2/3 [■*(•«»")„ - J} 


^vc(isr^) = S5 


2 A 


<2 \ 


1 


SAud COS ip ^ 3^ (■^^“)45 S A\i2l3 5 


SXp cos ifM = — 


— {sIm 

SV y \ 




3sv 
2 


<2 \ 


A"")™ - 


sxi = 


SX‘^P = 


SXu = 
5Xd = 
6 X 3 = 


2svu 

1 

2svd 

1 

1 


3?;^ sin 2/3 
t ! CaI \ cos/3 

■’A"").,- — 

^ ! Ca 2 \ sin/3 


<2 \ 


2v 




45 S 

; d\x 

\ / 45 


S(Mio) - —sin2/3(5(7W 




- cos^ /3S 


2A 


{«Ah”),,-“"’33'3A»")44} 




1 „ / 


44 


D.3 Leading two-loop effects 0{Y^\a.(^s)iY^f^) 

We follow [^[241: 


3Y‘ 

SXu= ‘ 


256A 


«4=“‘‘ 


256A 


16ff3 + ~ 3 3 cos^ I3Y^^ 




In - y 

mi 


+ [3 cos^ PY^^ + (3 cos" /3 + 1)^,"] (^In" ^ “ In" | 


1653 - o 5 " + 3 sin" /3y*" - 3 cos" /SF," 


Q" 


In — 5 - 

mi 


+ [3 sin" /3A," + (3 sin" /3 + 1)F*"] ( In" ^ - In" ^ 


m? 


mi 


Note that these leading effects are conveyed by the SM-fermion and gauge sector, as well as the doublet Higgs 
sector, so that the new-physics phases do not intervene. The contribution of sfermions or gauginos is merely 
reduced to the cutoff Q in the logarithms. 


D.4 Pole corrections 

Here we compute the shifts in the Higgs self-energies, which then allow to evaluate the pole corrections to the DR 
Higgs masses. We use the notation Af{p^) = f(jP) — /(O) for any function / of the external momentum. We are 
still working in the Feynmann gauge. 
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i) Contributions from SM fermions 

AHsosoip^) {N,Y^ + X^^X^JABFF{p,mu,mu) - 2ml[Xtxf^ - Xix/jABo(p,m„,m„) 

+ NCYj {{Xf^Xf^ + XI^X]^)ABff{p, md, m^) - - Xl^X]^)ABo{p, m^, m^)] 

+Y^ + X^aX^^,)ABMP, me, me) - 2ml{X^dXfd - X/,A/JABo(p, me, me)] } 

ABh+h-{p^) [{Y^ cos^ P + Y^ siv? j3)ABFF{p, mu, md) - 2YuYdmumdS\n2l3 ABo{p,mu,md)] 


+Y^ sin^ /? ABff{p, 0, me)] 


a) Contributions from gauginos and higgsinos 


AH so so {p'^) = 


IGtt^ 


(,'2 32 

— ABff{p, Mi,p) + —ABff{p, M 2 , p) + \^ABff{p, ms, p) 

X {XPuX^u + XfdXf^ + XluX]u + xl^x]^) 

+ 2 \)?ABff{p, P, p) + k^ABff{p, ms, ms)] (X^^Xf^ + 


1 


± 0'^ 3o^ 

Anff+R-(p^) =j^ I + —ABff{p,M2,p) + \^ABFF{p,ms,p) 


in) Contributions exclusively from the electroweak gauge sector 
1 7 


An5.ogo(p ) =y^^2 I 9 XIfyAB q{p, Mw, Mw) + 


,^2 


MzABo{p, Mz, Mz) 


X sin^ PX^X^u + cos^ PX^^Xf^ + sin P cos + X^^Xfu). 


AHh+h-{p^) =0 


iv) Contributions from the gauge / Higgs diagrams 
1 


ATtgOgoip ) = 


327r2 


/2 _|_ 2 

g^ABsvip, Mw, Mw) + ^ ABsv{p, Mz, Mz) 


+ 


X [sin" pxPuXfu + cos" + sin /? cos /?(X«Xf, + )] 

g'^ABsvip, Mw, Mw) [sin" /3A/„X/„ + cos" l3X)aX)a - sin/3cos/3(A/„A/d + 

+ g^ABsviP, mH±, Mw) [cos" /?(X« A/„ + + sin" l3{xP^Xf^ + 

+ sin p cos /?(A« A« + Xf^Xfu - A/„- X/^X/j] 


+ -—ABsv {p, mgo , Mz) x 


( ■!/■/ V "^ I V -^ V "^ 

^iu^ku — -^iu^ku + -^id^kd ~ -^id^kd 


AHu+H-iP^) =7 


) (x/„xf, - xfuXL + xfdXid - x/dXfd) I 

/ /2 _ 2\2 2 5 'I 

327r2 I g '2 g 2 ^^sv(p,m„±,Q) + ■X-,^J^ABsv(p,mH±, Mz) + ^ ABsv(p, mgo, Mvv) | 

x [cos" /?(X|L" + xL") + sin" /3(Xf/ + xU) - sin/? cos /?(X|Lxfd - xLxL)] 


2<7'"5" 


wj Contributions exclusively from the Higgs / Goldstone seetor 


Art ( 2\ 1 / SfH+H- S ? ir + ir “ a n i ' 7 I ^ S ? s “ s “ S ? S ® S “ AD /' ^ 

AngOsoip )=Y^\ 5 ‘ ™ "3 J " ™ABo(p,m^±,m^±) + - g ^ ^ ^g j r, r,^ABo(p,ms^,mgo) 


An 


H+H 


n=l,2 

/ 2 n 1 S°H+Hr S°H+H- 

(p)=Y^1^9- "fl™" 


ABo{p,mgo ,m^±) 


Tn=l,6 
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The tree-level Higgs-to-Higgs couplings are given in appendix |B.2.5| 
vi) Contributions from the sfermions 


^nsoso(P^) = 




E 


771,71 = 1 


' ABo ) + NcQ 


ABo{p, mf, ) 


S'lU'tUn S°-U*U„ 

1 9 

, S°E* E„ S'^E'‘EmAr, / , S^iV'JV S9]V*iV AD r ^l 

_l_g . m +g ' g ^ ABo{p,mj^,mj^)'> 


An 


H+H 


' (P^) = 


IGtt^ 


E 


H+U^D-n, H-D*U„ 


ABoip. 




,)+E5 


H+N*Em H-E*N 


ABoip. 


, "t'Err 


771,71 = 1 


The Higgs sfermion couplings are given in appendix |B.2.3[ 

E Simplified effective potential 

In this appendix, we study the simplified effective Higgs potential of Eq |33[ or more precisely the following and 
slightly modified version: 


=M^\S\^ + ^ + h.c] + Vo(|5|2) 

+ {Ml + Xf,\S\^)\H^\^ + (Mi + Xj.\S\^)\Hd\^ + + + h.c.] 


( 66 ) 


A3|H„r|HdP + A4|H„-ihdP 


■ Hdf + (A 6 e*^«|ih „|2 + + h.c. 

2 s 


This simplified potential is meant as an expansion of the effective potential - see Eq |29| - up to quartic order 
in the doublet fields. It slightly differs from Eqj^ in that the Zs-symmetry has been explicitly restored in the 
terms of the last line. Note that this way of restoring the Za-symmetry is just the simplest educated guess, while 
any additional factor /(|S'p, 5"*^) could intervene. Therefore, the factors of S/s appearing in the last line are 

just chosen as such because they will provide improved results numerically. Formally however, the associated 
corrections will remain of subleading order in the expansion in the doublet v.e.v.’s. 


E.l Matching the tree-level Higgs potential 

The tree-level Higgs potential (Eqj^ matches straightforwardly on Eqj^ 


M| = m| 

Ase^‘^^s = KA^e'-'^^ 

Vo(|5|^) = «:^|E 


Ml = m2 


Ml — mjj^ 


\u _ \2 _ \d 

Ap — A — Ap 

Aude^'^^-^ = XAxC'^^ 
X^ = 


K = = Ad 


A 3 — 


- 9 '^ + 9^ 


= 0 = = XrC'^^ 


E.2 Minimization conditions 

Ml = [{Aud cos + Ap S cos Pm) s + 3A6 cos pevl + X 7 cos pr^l] — - [Xps'^ + Xuvl + (A 3 -I- A 4 -I- A 5 cos ^> 5 ) 7 ;^] 

'^U 

Ml = [(Aud cos pA^d + Ap s cos pai) s + Ag cos p^vl + SAy cos prvl] — - [ApS^ -h Xdvl + (A3 -h A4 -h Ag cos P5)vl] 

Vd 

M/ = — (Vg) — A 5 S COS pas ~ Apu^ — Xpvl i^Aud cos pA.ud ff 2A ^s cos pm) —~— 

Aud sin pA^d = -Ap s sin pM + - [Ag sin p5VuVd - Ag sin pqvI - A 7 sin P 7 vl] 

Hssm(^^s = [AudSmpAud “ 2Ap ssinpAi) cs -3Ap sin^^M^— 
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O(w^) terms have been neglected. (Vg) = 45 f||5|2 2' 

E.3 Higgs mass matrices 

Charged Higgs in the base {H^, ): 

(At^±)ll= 0 ±— {■^H±}22 ^ 

= Re cos * cos (Pm) s — (A4 + A5 cos (^5) VuVd + Ag cos (fieVu + ^7 cos (frVd ( 67 ) 

Neutral Higgs in the base (h°, a°, a°, a°); 0 {v*) terms are neglected.: 

= [{AudCosp>A^^ + XpS cos p>m) s - SAecosvJewS + X-j cos ‘p'jvf\ — + 2 A„w^ 

(Mpo)^^ =- [{Aud cos (pA^a + XpS cos (Pm) s + 3 Xe cos (pevl +3X7 cos iprvj] + 2 [A3 + A4 + A5 cos vJs] 

(Mpo )22 = [{Aud cos (PA^^ +XpS cos Pm) s + Ag cos pqvI - 3X7 cos P7vl] — + 2Adz;^ 

'^d 

= - Aud cos pA^a + 2 Xp SCOSPM + ^ ( 3 Ag COS(/?gw 2 + A7 COS - 2A5 COSp^VuVd) Vd + 2 ApSW„ 

{Mjjo )23 = - Aud COS pAui + 2 Ap S cos Pm + ^ (Ag COS Pevl + 3X7 cos P7vl - 2A5 COS P 5 VuVd) Vu + 2 XpSVd 
{Mpo )33 =S 24 s cos PAs + 2 s^ (Vg ) + Aud cos PAuA 

(■^770)44=^07“ ! (-^770)45 =Hg ; (7Wpo)55 = Hg — 

^d 

Hg = {Aud cos Pa^u + Xp S cos Pm) S - 2A5 cos p^VuVd + Ag cos PqvI + X7 cos P7vl 
= Aud cos pAud - 2 Ap S cos Pm + ^ (Ag cos Pevl + X7 cos P7vj - 2X5 cos P5VuVd) Vd 

{Mjjo)^g = Aud cos P Aud - 2 Xp SCOSpM + ^ {XeCOSpevl + X7COSp7vj - 2A5 cos Vu 

(Atpo)gg = - SsHs cos (^As + {AudCospAud + 4 :XpS cos Pm) 

(Atpo)44 =( 2 Ag sim^gUt, - Ag sinp5Vd)vd 
(Atpo)45 = ( 2 Ag sinV 5 g 7 ;„ - Ag sinV 3 gt!d) Vu 

{Mpo)^Q = - 3 Xp s sin Pm + ^{ 2 Xe sin Pevl- X5 sin p^VuVd) Vd 

(Mpo )24 = (2A7 sin p7Vd - Ag sin p^Vu) Vd 
{Mpo )25 = (2A7 sin p7Vd - Ag sin p^Vu) 

(Atpo)2g = - 3 XpSsinpM+ ^{2X7sinp7v‘^-X^sinp^VuVd) 

=Vd I^Ap ssmv 5 M-sm(/?g 

{Mpo /gg =Vu I^Ap s sm PM - — Ag sm p^ 

(Atpo)3g = 4 Ap sinpMVuVd 
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E.4 Trilinear Higgs couplings 


qO c^O oO 1 


^u'^u 


V / t,J,K 


' ' '^■>3,1^ 
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+ A5 cos (P5 
+ A5 sill (^5 
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(n»)"" + (n-')'-") + va + (d 

:n'’):j:'-(n"):,;“- 2 (n-)t-r)+.„((i 

n')r;‘:'-(n'')“"-2(n'')?:rn+»A(i 


i,j,k 

+ 


+ Ae sin ipe 
— Ay COS (fiT 
+ Ay sin(^y 


+ A^ 


.((n")Sr 2 m:;?) 

u,, - - 2 (n-)*".-') +»'(<"')S:: - (n'-)S;: - 2 (n^-^)] 

•V r f ^-r-r‘^\d,u,u /^A\d,u,u ^ /^A\U,u,d\ / /^q\U,u,u /^a\U,u,u 

A6cosw[»..(3(n®).j,j +(n'*)y,^ -2(n^),j„)+»j((n''),.^ +(n-").j,t 

' - ^»4^<"")J»+2(nn:«+m:rr)+»4<"')“rr+(n')"rr; 

((n^)*" + (n-*)"-) +.. (3 (n-)”" + (n-^)"" - 2 (n-)".”) 

. ((nni" + (n-)Z::) + (3 (n'-)“" + 2 (nO": + (n')"") 

_4..cos^. [(n^)-;-:' - (n^)--- _ (n^)-" _ (n^)*"- 

A? cos^M [ 2 , ((n^):-:- - (n-)'*:* + (n-)';-- + (n--)J~) 

++2 (n%:: - ((n«” 

(3 (n')'*^ - 3 (n'-)-;-- + (o'-)-;- + 

+[> ((n*):;"r+(n^)":)++(n^):--) 

+ [. ((n*)-;" + (nAj-’i") + „ ((nA)?;; + (nA) J;“) 

+ [(nA)-;" - 3 (nA)'"l + ?,3 ^ 


"*)“: +2 (nA)“;-(nA)“--)] 

iunt:’:)-"-f{(ur:z-Hnn:x 

«3((n0::;:: + 2 (n'>):;;;:-(n'-)";")] 


S"HJH 7 
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XuVuX^X^^X^^ + XdVdX^X^^X^^ + As VuX^X^^X^^ + Vd cos^ l3X^X':t„X, 


rC 

•-ju^ 


C ■ 

ku_ 


+ 2 [~(a ^4 + A 5 cos (ps) [vuX^ + VdXj^) + As sin (p^ {vuXl^ + VdX^)] + X^^X^^) 

+ [Ae cos ipeVuXj^ + Ay cos ip^VdX^] [X^^X^^ + Xj^AT^) 

+ [-As cos (/?6 {vuX,^ + VdX,^) + Xe sin ipQ {vuX^^ + VdX^)] Xf^X^^ 

+ [-AyC 0 S(/ 9 y {vuX^ + VdX^) + Ay sin^Sy {vuXld + VdXQ] Xf^XjXd 

+ [XpXgxZ + XU^dXZ] sX^s 

+ 2 [(^^“d COS (fA^d + ‘^Xp s COS (Pm) X^ + 3Xp s sin (pMX-g\ (X^X^ + Xj^X^) 

+ ^ [(A4 - As cos (/9s) {vuX^d + VdXk) - As sin (/9s {vuX^ + VdX^) 

+ 2 (As sin ipQVuX^ + Ay sin iprVdX^) 

+ {Aud cos ipA^d - 2A7s cos pm) X^ + X^s sin pmX^] (X^Xf^ - X^X^) | 
We omit the quartic couplings. 


F Extension to Zs-violating terms 

A departure from the Zs-conserving NMSSM is motivated by cosmological considerations (Domain-Wall problem). 
Turning to the most general singlet extension of the MSSM, the solution to the /i-problem becomes less obvious, 
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as several dimensionful parameters now enter the superpotential. However, this difficulty is lifted when these 
supersymmetric masses appear as a by-product of the supersymmetry-breaking mechanism - see e.g. 10 11 . This 
motivates the inclusion of Za-violating terms. Here we discuss how our results can be extended to this more general 
version of NMSSM. 

Eqs(^and[^are supplemented with the following operators: 


AIT = ne'^'^Hu ■ Hd + ^Fe"'^^S+ 

/ 2 

- ■ Ha + S + + h.c. 


( 68 ) 

(69) 


However, without loss of generality, the /i-term (for instance) can be set to 0 by means of a shift of the superheld 
S and a re-dehnition of the other dimensionful parameters: 


A 




(70) 


•<— ^'e 


L p'-'fr, 


mo e 


A 

2 / 
A 


^ _ 2 ^\^q^(V 2 +<I>u-Vx) 

2 A 


A 2 

,2 


A^ A A 

This choice simplifies significantly the corrections to the Za-conserving case and we thus use this freedom in the 
rest of this appendix: /r = 0. We then have five new complex parameters with respect to the Za-conserving 
NMSSM. Note that these only affect the singlino and the Higgs sector so that most of what we derived in the 
context of the Za-conserving NMSSM remains valid. 

The singlino mass-entry in the neutralino mass-matrix - Eq{T^- is changed to -I- 2ft:e®‘^'‘s and this is 

the only modification in the gaugino / higgsino sector at the order of our calculation. Concerning the neutralino 
loop corrections to the Higgs sector, our results of appendix |D. 2.4 can be extended to the Za-violating case by the 
simple substitutions: 2«;e*‘^''5' —)■ -I- 2ft:e®‘^'"S' and —)■ + 4/i'«:Re -I- iK^\S\'^. 

The modihcations in the Higgs sector are more substantial. The tree-level Higgs potential of Eqj^ receives the 
following additions: 


AVh = 


-k A -k {H+R- - H°H^) + h.c. 




e 




S- 




S'^+ K^i'e^'^^'^-'^F■)S\S\^ + h.c. 

+ eF + As? (71) 

Correspondingly, the Higgs mass-matrix elements of Eqj^and ^ are modified as: 

AH = m\ cos y >3 -I- \^F cos((pA — ^f) + As/i^ cos(:^a — 

a(m^o ),2 = -ah 
a(m^o )43 = ah 


Am?,, = AH^i±^ 


VuVd 
Vd 


(72) 


A{M]jo),, = An^ 


HZ 


A{Mlo),, = An^ 


A (Mho)^^ - AH — 


Vd 


A = -XvdfJ.' cos{(px - 


A (A4^o) 23 = -Awu/r' cos(v?A - ^Pf') 

VuVd 


A(A4fo)33 = -- [^s COS if s + ip h' cos{ipy,> - pf)] + h' 3«;scos((^«, - py,:) + \— —cos((/?a - Pf') 


A = -Xvdfi' cosipx - Pf,>) 


A (Mho)^q = -Xvuh' cos{px - Pf') 


A(Mho)qq = -- l^s cosps + cos{pf,< - v5f)] + p' -kscos{pi^ - v?^/) -k cos((/?a - Pf') 


- 2 


m'g cos Pm' + 2 k^f cos((^k — Pf) 
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A , 15 , 24 , 25 ® 

A = - [mgsintpa + A^f sin((/2A - V^f)] y A = - [mgsintpg + A^f siii((/2A - V^f)] ^ 

A (Mjjo)-^Q = - [mgsint/jg + 2As/i'sin((^A - + A^f sin(95A - V^f)] y 

A (Mjjo)^Q = - [mgsintpg + 2Xsn'sm{(px - + A^f sin((pA - Vf)] ^ 

2 2 

A {Mho)^q = - [Cssintps + sin{(pf_,, - ipp)] + m'g sinipm' + 2K^Fsin{ip^ - (pf) 

+ [™3 sinyjg + A^Fsin((^A - ^f) + 2Xs fi' sm{px - <ftj,')] 

Finally, the trilinear Higgs couplings of appendix |B. 2. 5| receive the following changes: 


s"h+h: 


=;^{ 


Xp' cos((/?A - + { -[^3 sin ‘Ps + A^f sin((pA - Pf)\ + 2A/i' sin((pA - P^i')] X 


— % 


X (Af„A, 

Xp! cos{p\ - + - [ml sin(/jg + A^f sin((/?A - Pf)] Xl 


c 

kd 


■xf,xg) 


(73) 


f \ 

K^ju^kd ~ ^jd^ku) 


gS°S°S° =;^| - xg cos{px - Pf,>) [nf^r'^ - + Kg cos{p^ - p^') [n 


S sss 
ijk 


n 


A sss] 

ijk J 


- {ml sin p 3 + X^F sin(i^A - Pf)) 


P sud 


- 2Xg sin((^A - pA 


H: 


P sud 
ijk 


+ 


3s2 


n 


n: 


P uds 


ijk ' ^^ijk 
'^u'^d ri-r/ sss " 


ttP dus tt/ sud 

^^ijk ^^ijk 


VuVd 

3s2 


[3n^r* 


TTi SS^ 

~ ^hjk 


- ^ijk 


*] - o^g sin((/j,, - pA^ifg 


- ^ Q[^ssinv?s + ^FAi'sin(v 2 ^/ - (^f)] + sin [nf/jf" - 


Up to these modifications, the calculation of Higgs loop corrections to the Higgs masses - see Appendix |D.2.5 
{in, iv) - remains valid. On the other hand, reconstructing the contributions to the effective potential ~ i.e. 
the effective couplings - in the same fashion as in the Zg-conserving case - see Appendix D.2.5 {v) - becomes 
problematic: indeed, no argument of symmetry allows to reduce the number of terms in the generic potential 
for two Higgs doublets and a singlet and the latter contains too large a number of parameters in view of the 16 
independent mass entries. Refer to 28 for a discussion concerning the generic Higgs potential. It is possible to 
use an expansion of the Higgs mass matrices in terms of the doublet fields, however. 
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